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Abstract 

We refine and develop the inverse scattering theory on a lattice in 
such a way that the Ablowitz-Ladik lattice and derivative NLS lattices 
as well as their matrix analogs can be solved in a unified way. The 
inverse scattering method for the (matrix analog of the) Ablowitz- 
Ladik lattice is simplified to the same level as that for the continuous 
NLS system. Using the linear eigenfunctions of the Lax pair for the 
Ablowitz-Ladik lattice, we can construct solutions of the derivative 
NLS lattices such as the discrete Gerdjikov-Ivanov (also known as 
Ablowitz-Ramani-Segur) system and the discrete Kaup-Newell sys- 
tem. Thus, explicit solutions such as the multisoliton solutions for 
these systems can be obtained by solving linear summation equations 
of the Gel'fand-Levitan-Marchenko type. The derivation of the dis- 
crete Kaup-Newell system from the Ablowitz-Ladik lattice is based 
on a new method that allows us to generate new integrable systems 
from known systems in a systematic manner. In an appendix, we de- 
scribe the reduction of the matrix Ablowitz-Ladik lattice to a vector 
analog of the modified Volterra lattice from the point of view of the 
inverse scattering method. 
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1 Introduction 



The cubic nonlinear Schrodinger (NLS) equation [HE] is probably the most 
prominent example of an integrable partial differential equation in 1 + 1 
space-time dimensions. The inverse scattering method for the NLS equation 
devised by Zakharov and Shabat [TJ was reformulated by Ablowitz, Kaup, 
Newell and Segur [2111] in a user- friendly and broadly-applicable manner. 
Since then various extensions of the NLS equation have been obtained within 
the framework of the inverse scattering method. Among them, we mention 
two kinds of extensions: (i) space-discrete NLS system^] wherein the spatial 
variable is discretized JBJ[T] and (ii) derivative NLS systems wherein the non- 
linear terms involve differentiation with respect to the spatial variable [5HT2]. 
For these extensions, the inverse scattering method can be applied on a case- 
by-case basis, but its application requires more steps and is apparently more 
complicated than that for the original NLS system [3111] ■ 

The main objective of this paper is twofold. First, we develop the inverse 
scattering method on a lattice and correct the widespread impression that it is 
essentially more complicated than the inverse scattering method on the line. 
Second, we show that the derivative NLS systems can be solved using the 
inverse scattering method for the NLS system; however, this paper focuses on 
the space-discrete caseH To be specific, we consider (a matrix generalization 
of) the Ablowitz-Ladik lattice [B] that is an integrable space discretization 
of the NLS system. The inverse scattering method for the Ablowitz-Ladik 
lattice reported in the existing literature involves some onerous processes pe- 
culiar to the discrete case; in fact, they are redundant and can be avoided. 
For this purpose, we only need to start with an eigenvalue problem that is 
trivially equivalent to the Ablowitz-Ladik eigenvalue problem up to a simi- 
larity transformation and inversion of the spatial coordinate. Then, all the 
key quantities such as the scattering data become even functions of the con- 
ventional spectral parameter; thus, we can use its square as a new (and more 
essential) parameter. This considerably simplifies the subsequent computa- 
tions. Moreover, by the inversion of the spatial coordinate, we no longer 
need to normalize the "integral kernels" of the linear eigenfunctions (Jost so- 
lutions) to express the potentials in the Ablowitz-Ladik eigenvalue problem 
explicitly. This is in contrast to the conventional approach wherein one has 
to first introduce the "integral kernels" and then normalize them; in the lit- 
erature (see, e.g., [TB1[T6]). they are usually denoted as K(n, m), K(n, m) and 
K(n,m), K(n,m), respectively. Thus, we can successfully refine the inverse 

lr The problem of how to discretize the continuous time variable in integrable space- 
discrete systems has a long history, see [5] and references therein. 

2 Relevant results on continuous derivative NLS systems can be found in [LH[T^] . 
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scattering method associated with the (matrix) Ablowitz-Ladik eigenvalue 
problem; both the potentials and the linear eigenfunctions are determined 
from the scattering data through a set of linear summation equations in the 
most transparent manner. 

In our previous paper [UJ , we proposed a systematic method of generating 
new integrable systems from known systems through inverse Miura mapsjj 
As a result, two derivative NLS systems, namely, the Gerdjikov-Ivanov (also 
known as Ablowitz-Ramani-Segur) system pUHH] and the Chen-Lee-Liu 
system [9], were constructed from the Lax representationE] for the NLS sys- 
tem; the same prescription applies to the space-discrete case. Thus, the 
inverse scattering method for the Ablowitz-Ladik lattice can also provide 
the solutions of the space-discrete Gerdjikov-Ivanov system and the space- 
discrete Chen-Lee-Liu system in a unified way. In this paper, we propose yet 
another method of generating new integrable systems from known systems. 
In particular, by applying this new method to the Ablowitz-Ladik lattice, 
we obtain a lattice system that is essentially equivalent to the space-discrete 
Kaup-Newell system studied in [SUIET]; its solutions can be expressed in 
terms of the linear eigenfunctions associated with the Ablowitz-Ladik lat- 
tice. Thus, the space-discrete Kaup-Newell system can also be solved us- 
ing the inverse scattering method for the Ablowitz-Ladik lattice; note that 
among the derivative NLS systems, the Kaup-Newell system [8] is the most 
important for physical applications. 

The main body of this paper is organized as follows. In section 2, we 
start with the Lax representation for the Ablowitz-Ladik lattice; using its 
linear eigenfunctions, we derive two derivative NLS lattices. First, we apply 
the method proposed in [UJ and derive the space-discrete Gerdjikov-Ivanov 
system. Second, we propose a new systematic method and obtain the space- 
discrete Kaup-Newell system. We can also derive and solve the space-discrete 
Chen-Lee-Liu system, but we do not present it in this paper; the interested 
reader is referred to appendix B of [17] and section 3 of [20]. In section 3, 
we present a streamlined version of the inverse scattering method for the 
Ablowitz-Ladik lattice. In section 4, we combine the results of sections 2 
and 3 and show that the derivative NLS lattices can be solved by the inverse 
scattering method associated with the Ablowitz-Ladik eigenvalue problem. 
Their multisoliton solutions can be derived from the linear summation equa- 
tions in a straightforward manner. The last section, section 5, is devoted to 
concluding remarks. 

3 The original Miura map transforms the modified KdV equation (or, more generally, 
its one-parameter generalization called the Gardner equation) to the KdV equation [TS] . 
4 The term was coined after Lax's work on the KdV hierarchy |19j . 
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In this paper, we consider the general case where the dependent variables 
take their values in matrices [3122] in such a way that the operations such 
as addition and multiplication in the equations of motion make sense. In 
appendix [A] we consider the reduction of the matrix Ablowitz-Ladik lattice 
to a vector analog of the modified Volterra lattice and discuss the effect of 
the reduction on the scattering data; this considerably refines our previous 
results given in [22JI21]- 



2 Ablowitz— Ladik lattice and derivative NLS 
lattices 

2.1 Lax representation for the Ablowitz— Ladik lattice 

We start with the matrix Ablowitz-Ladik eigenvalue problem written in the 
nonstandard formjfl 



zl 

~ 1 R r , 



*l,n+l 
^2,n+l 



(2.1) 



where z is a constant spectral parameter and / is the identity matrix of 
arbitrary size. For simplicity, we assume that all the entries in ( 12. ip . such 
as the potentials Q n and R n , are I x I square matrices. It is also possible to 
consider the more general case where Q n is an l\ x l 2 matrix and R n is an 
l 2 x li matrix; however, the results for that case can be easily obtained by 
setting some rows and columns in Q n and R n as identically zero. 

The time evolution of the linear eigenfunction can be introduced in such 
a way that it is compatible with the eigenvalue problem (12. ip . The most 
illustrative example is given by 

-z 2 + l)bl + bQ n _ x Rn -z 2 bQ n ^ - aQ n 







. ^2,n _ 





— bR n — z aR n . 



1 



')al + aR n _iQ r 



(2.2) 

Here, a and b are arbitrary scalar constants; in fact, they may depend on 
the time variable t in an arbitrary manner (see, e.g., [251127]). but we do not 
discuss it in this paper. The compatibility condition for the overdetermined 
linear system, (12. ip and (12. 2p . with the isospectral condition z t = implies 
the time evolution equations for Q n and R n : 

Q n ,t - aQn+i + bQ n -i + (a - b)Q n + aQ n R n Q n+ i - bQ n _ x R n Q n = O, (2.3a) 

Rn,t — bR n+ \ + aRn-i + (b — a)R n + bR n Q n R n+ i — aR n _iQ n R n = O. (2.3b) 



5 Actually, in the simplest 2x2 matrix case, (|2.1I) can be identified with a special case 
of the eigenvalue problem studied in [33] . 



5 



We call (J23D the (matrix) Ablowitz-Ladik lattice/system gUE]; (EHJ) and 
(12. 2p comprise its Lax representation. The symbol O on the right-hand side 
of the equations implies that the dependent variables can take their values in 
matrices. In fact, there exist infinitely many ways to define such an isospec- 
tral time evolution (cf. [281 ES]), depending on the choice of the temporal 
Lax matrix in (12.21) as a Laurent polynomial in z 2 ; they provide the positive 
flows of the Ablowitz-Ladik hierarchy (cf. appendix A of J2D] for the neg- 
ative flows) and each of them is uniquely determined by its linear part or, 
equivalently, the dispersion relation. 

2.2 Space-discrete Gerdjikov— Ivanov system 

In this subsection, using the method proposed in [17], we derive the space- 
discrete Gerdjikov-Ivanov system from the Lax representation for the Ablowitz- 
Ladik lattice. The result is essentially the same as that given in [17J, but we 
restate it here for the self-containedness and readability of the paper. 

We consider a 21 x I matrix- valued solution to the pair of linear equations 
(12.11) and ( 12. 2 p such that \l/i in is an / x / invertible matrix. Then, in terms of 
the I x I matrix P n : = ^ 2 ,n^rn> (EU) and O can be rewritten as a pair of 
discrete and continuous matrix Riccati equations for P n , 

R n = fiP n - P n+1 + fiP n Q n P n+1 , (2.4a) 

P n>t = -bR n - [i^aRn-i + (1 - n~ l )aP n + l)bP n 

QnPn bP n Q n —iR n -\- [ibP n Q n —\P n -\- dP n Q n P nj (2.4b) 

where /z := z 2 . The first relation (I2.4ap defines the Miura map (Q n , P n ) H- (Q n , R n ). 
Using (I2.4ap . we can eliminate R n and -R n -i in (I2.3al) and (12.4bl) to obtain 
a closed system for (Q n ,P n ), i.e., the space-discrete Gerdjikov-Ivanov sys- 
tem [20]: 

' Qn,t - aQn+i + bQ n -! + {a - b)Q n + aQ n (fiP n - P n+1 ) Q n+X 

— bQ n ^i (flP n — P n +l) Qn + O'^QnPnQnPn+lQn+1 ~ bfiQ n -lP n Q n Pn+lQn = O, (2.5a) 

< 

P n ,t - bP n+1 + aP n _i + (b - a)P n - bP n (Q n -i - nQ n ) Pn+i 
k + aPn-i (Qn-i — ^Qn) Pn + b[iP n Q n _\P n Q n P n +\ — afj,P n -iQ n -iP n Q n P n = O. (2.5b) 

2.3 Space-discrete Kaup— Newell system 

In this subsection, we propose yet another method of generating new inte- 
grable systems from known systems using a fundamental set of linear eigen- 
functions. The method is applicable to a matrix Lax representation of ar- 
bitrary size, but for brevity we describe it in the simplest case of a 2 x 2 
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Lax representation as well as its block-matrix generalization involving two 
potentials. In this paper, we consider the space-discrete case (see [H] for the 
continuous case). 

Suppose that a lattice system admits the Lax representation 



1,71 



(?) 
2,71 



Ln,n 
£21,71 



L\2,n 
L22,n 



^ 1,71+1 



(?) 
2,71+1 



(2.6a) 



71 

(i) 

2,71 



M 2 i n 



Mi2,„ 

M 22 „ 



1,7J 



(7) 
2,7! 



(2.6b) 



where all the entries are assumed to be square matrices of the same size. 
Because the method requires a full set of linearly independent eigenfunctions, 
the superscript with j = 1 or 2 is used to designate the two eigenfunctions. 
We apply a gauge transformation defined using one eigenfunction to the other 
eigenfunction as 



(2) 



(2) 
2,n 



K+ 



o 

2,71 



(2) 



(2) 

2,71 



l,n l,n 



Then, the Lax representation (12. 6p is transformed to the degenerate form: 



J - W l,n L 12,ri*2,n+1 



(i) 



(1)-1 



12,71*2,71+1 



71+1 



^(1)"1^(2) 

V 1,71+1 ^ 1,71+1 

^(1)"1^(2) 

V 2,71+1 V 2,71+1 



(2) 



(i) 



(i)- 

2,7! 



1,71 1 



(2) 



2,n 



^(2) 

2,71 



Indeed, there are only 2 + 2 independent quantities in these Lax matrices. 
Thus, we can express them in terms of the components of the linear eigen- 
function as 



a)- 



-^12,71* 



2,71+1 



^2,71 "^21,71* 1,71+1 



and 



qrW-i-M vl/ (1) 

^)- l M H/ (1) 
^2,n M 21,n^i, n 



*l,n *l,n *l,n+l*l,n+lJ i *2,n+l*2 ) n+l * l,ra+l * l,n+l 

(2.7) 

^(1) -1^(2) _ ^(1) -1^(2) \ / ^(1) -W2) _ ^(1) -1^(2) 
^2,n w 2,n v 2,n+l * 2,n+l J l*l,n+l *l,n+l ^2,71+1^2,71+1 

(2.8) 



- • (2-1°) 



-1 
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Note that ( J2,6ap also implies the two important relations: 

*S"%i,»M5Ui = 1 - ( 2 - n ) 

V»*£Ui = I - (2.12) 
Typically, n and £22,™ are ultralocal functions of L 12 n and L21.ru such as 

L\l )n = Oil + /3L 12)n I/21,n; ^22,71 = 7-^ + 0-^21,71-^12,71) 

where a, /3, 7 and 5 are scalar functions of the spectral parameter. Thus, we 
can try to solve fl2~TT|) and fF2~T2l to express ^"^K+i and ^"^S+i 
in terms of ^^-k^n^n+i 

and ^l/^n 1 -^2i,n^'i^+i- This is relatively easy 
if either Ln.„ or L22,n is a constant scalar matrix, e.g., (3 or 5 vanishes in the 
above example. Then, using (12 .7p and (12.81) . we can also express 



*l,n -^12,71^2,71 - W l,rt ^12,n V 2,n+l ^ W 2,n W 2,rt+1 

^1,71+1^12,71^2,71+1 - I W l,n W l,n+1 J V l,n ^12,n V 2,n+1> 
^(l)-lr vl/(l) _ *>(1) -1 r ^(1) /^(l)-l^(l) \ _1 

y/ 2,n+l- L 21,n^l,„+l - [^2,n V 2,n+1 J ^2,71 ^21,n*l,n+l> etC - 

recursively in terms of 

:= ^ : =^"^S- (2-13) 

In general, M 12i „ and M 2 i in are local (but not ultralocal) functions of Li2, n 
and L2i, n - Therefore, with the aid of the above relations, (12.91) and (I2.10p 
can be rewritten as a closed lattice system for u n and v n . 

Let us illustrate the method using the matrix Ablowitz-Ladik lattice (12. 3p 
as an example. By setting 

-^ll.rt = Zl, L 2 2,n = z /, 

(12. lip and (12.121) provide the useful relations 

*£+!*£ = z{l~ ^VL^Uiy 1 , (2-14) 



V 2,n+1^2,n — z 



(/ - *^%i,n*£+l) • (2-15) 
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The Lax representation, (12.1 1) and (12.21) . implies simple relations between 
off-diagonal elements of the temporal and spatial Lax matrices, i.e. 

M 12j n + z^aL^n + zbL 12: n-i = O, 
M 2 i, n + zbL 21>n + z~ x aL 21 ^ x = O, 

which can be rewritten as 

\&W _1 flyf Vl/ (1) 4. r-^n^^T H> (1) fl/W "W 1 ) 

v ln ik/i 2in w 2n + z aw ln Li 2)U m 2n+l m 2n+1 m 2n 

+ ^ 0W l,n y/ l,n-l W l,n-l- L 12,n-1^2,n ~ U > 

n> (1) 4 _1 r H/ (1) v&W-^W 

+ 2 «^2,n y/ 2,n-i y/ 2,n-l- L 21,n-l^l, n — V . 

Substituting (gUD and (r2~T5j) and subsequently (T2Tj) - ([2TTUj) into the above 
relations, we arrive at a closed system for ^i^ -1 ^ „ and ^f 2n ~ 1 ^ r 2n- 



Proposition 2.1. Consider two linearly independent solutions of (12.11) and 
( 12.21) and write them as in (12. 6p with j — 1 or 2. Then, u n and v n defined in 
(I2.13P satisfy the following system: 

Un,t + ~ (U n ~ U n+ i) [I + (V n - U n )~ l {u n - U n+l )] 1 

+ fib [I- (tt n _i -u n )(v n -M„,) _1 ] (Un-i ~ U n ) = O, (2.16a) 

< 

fn,t + (fn - Un+l) [/ + ( M n ~ Vn)" 1 ^ - V n+1 )] 

+ - [I - [v n -i - v n ){u n - VnY 1 ] 1 (v n -i - f„) = O, (2.16b) 
V /i 

where /i = z 2 . 

Remarks: 

(i) The system (12.161) with fib = (a/fi)* allows both the complex conjuga- 
tion reduction v n = au* n and the Hermitian conjugation reduction v n = au^, 
where a = ±1. In addition, by setting v n = —u n , (12.161) with fib = a/ 1 \i re- 
duces to a single matrix equation, 

Un,t = U n [{U n+1 + U^ 1 - (u n + U„_l)~ ] U n , 

up to a rescaling of t. In the scalar case, this belongs to Yamilov's list of 
Volterra-type lattices in [30]; in the matrix case, it allows further reductions 



9 



to various multicomponent systems (cf. [3T]). 



(ii) The system (I2.16P provides a space-discrete analog of the system studied 
by Svinolupov and Sokolov [32j ; their system gives a matrix generalization of 
the Heisenberg ferromagnet model written in a two-component form [33|l3^] . 
Indeed, ( 12.1 6p in the scalar case is closely related to the lattice Heisenberg 
ferromagnet model [35] and its simplest higher symmetry [3U|I3"T]. 

(iii) The system (I2.16P in the scalar case appeared in the recent paper [38] 
(also see [39] ). In this context, it is natural to rewrite (I2.16P as a two- 
component system for the pair of variables u n and v~ l (cf. (4.8) in [33]). 

In (I2.16p . the two variables u n and v n interact with each other through 
the quantity (v n — M n ) _1 , which can be used as a new dependent variable. 
Indeed, a direct calculation shows the following two propositions. 



Proposition 2.2. Let u n and v n satisfy the system (12.161) . Then, the new 
pair of variables q n and r n , 



q n := A+u n (= u n+ i - u 



n) i '"n • \Vn ^n) i 



satisfies the space-discrete Kaup-Newell system [20J: 



{I-q n r n ) 1 q n + fxb (I + q n -ir n ) 1 Qn-i 



where denotes the forward difference operator. 



fib (I + r n g n _i) 1 r n + - (/ - r„_ig„_i) 1 r n _ x 



O, (2.17a) 
O, (2.17b) 



Proposition 2.3. Let u n and v n satisfy the system (12.161) . Then, the new 
pair of variables q n and r n , 

q n ■= (v n ~ u n )~ l , r n := A+w n „i (= v n - v n - X ) , 

also satisfies the space-discrete Kaup-Newell system (12.171) for q n and r n . 
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3 Inverse scattering method for the Ablowitz— 
Ladik lattice 



3.1 Revisiting the Ablowitz— Ladik eigenvalue problem 



In this section, we describe the inverse scattering method associated with the 
matrix Ablowitz-Ladik eigenvalue problem (12. ip . 



zl 

~ l Rr, 



ZQ T 



(3.1; 



Here, the potentials Q n and R n are assumed to decay sufficiently rapidly at 
spatial infinity: 

(3.2) 



lim Q n 

n—>±oo 



lim R n = O. 

n— >±oo 



The matrix generalization of the Ablowitz-Ladik lattice [6] first considered 
in the early 1980s [7J is still a topic of interest in discrete integrable sys- 
tems (see [20JI101I1T] and references therein). In our previous papers [23|[2¥j. 
we presented the inverse scattering method for the matrix Ablowitz-Ladik 
lattice while assuming some symmetry conditions on the potentials Q n and 
R n . Here, we remove such assumptions and consider the general case of I x I 
square matrices Q n and R n ; recall that the results on rectangular matrix po- 
tentials can be obtained by setting some rows/columns in Q n and R n as zero. 
The inverse scattering method reported here bypasses some redundant com- 
putation and consideration contained in the existing literature on the same 
subject, so we believe that this is the most streamlined version. The results 
in the previous work [23l|23] can be reproduced by imposing some reduction 
conditions on the scattering data; this is briefly sketched in appendix [A] In 
addition, we fix some minor inconsistencies in [231121], though they do not 
affect the main results of these papers. 

All the flows of the matrix Ablowitz-Ladik hierarchy are associated with 
the same eigenvalue problem (13. ip . so they can be solved together by the 
inverse scattering method. However, in the following, we concentrate on the 
matrix Ablowitz-Ladik lattice (12 .3p to illustrate the method in an easy-to- 
read manner. We stress that in contrast to other methods of obtaining special 
solutions, the inverse scattering method can provide the general solution 
formulas. Moreover, the method can determine not only the potentials Q n 
and R n but also a fundamental set of linear eigenfunctions, which will be 
used in section 4. 
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3.2 Jost solutions and relevant quantities 

To analyze the general case of the matrix potentials Q n and R n , we consider 
the adjoint equationjf] 



[ $i, n+ l $ 2 ,n+l ] = [ $l, n $2,n ] 

Indeed, a discrete analog of Lagrange's identity. 

[$l,n $2,n] 

[ $l,n+l $2,n+l ] - [ $l,n $2,n ] 



(3.3) 



^2,n 



Z.T zQ n 



^2,n+l 

zl zQ„ 



} 









implies that (13.1 p and (13. 3p can be said to be adjoint to each other. Thus, 
we can introduce an I x I matrix function W[- , •] for a pair of solutions to 
(HUD and (|33|| as 



#2,n 



$l,n*l,n + $2,^2, 



which is n- independent: 

W[Q n ,* n ] = W[$^ 1 ,* r *. 1 ). 

In addition to the rapidly decaying boundary conditions (13. 2p . we assume 
that the spatial Lax matrix defining the eigenvalue problem is invertible: 



det (I - Q n R n ) ^ 0, VnG Z. 



(3.4) 



Because log [det (/ — Q n R n )] is a conserved density for the matrix Ablowitz- 
Ladik hierarchy [Tj, this assumption is preserved under the time evolution. 
Thus, a set of column- vector (or row- vector) solutions to the eigenvalue prob- 
lem (13. ip (or (13.31) ) that are linearly independent at some lattice site, say 
n = no, remain independent for all n 6 Z. 



6 Note that if we consider a square matrix solution to the eigenvalue prob- 
lem \&„ = L n (z)^ n +x, then its inverse <I> n := satisfies the eigenvalue problem 
$„+i = $ n L n (z). 
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We introduce Jost solutions 



[z) and ip n (z), ip n (z) at a fixed time 



that satisfy (13. ip and the boundary conditions, 



I 

o 



and 



z n i) n -> 



> as n — >• — oo 



(3.5a) 



I 

O 



> as n — > +00. 



(3.5b) 



The time evolution of the Jost solutions will be considered in subsection 13.41 
Note that the overbar does not mean complex conjugation in this paper. Sim- 
ilarly, we introduce adjoint Jost solutions <^ d (z), <j>^(z) and ^ d (z), ip^{z) 
that satisfy (I3.3P and the boundary conditions, 

z n^d _^ J q _j j 
z -n^d — [JO] 

and 

z -n^,ad — ^ [ J O ] 



as n — >• 



-00 



(3.6a) 



n„/,ad 



->■ 



[<"] 



as n — > +00. 



(3.6b) 



Because the I + I (= 2/) columns of the Jost solutions ip n and form a 
fundamental set of solutions to the eigenvalue problem (13. ip . we can set on 



the unit circle \z\ = 1 as 



<p n (z) = 1p n (z)A + 1p n {z)B, 
<j> n (z) = i) n {z)B - 1p n {z)A. 



(3.7a) 
(3.7b) 



Here, A, B, B and A are n-independent I x I matrices, which depend on 
the spectral parameter z and are called scattering data. According to the 
asymptotic behaviors of the Jost solutions (I3.5p - fl3.6p . we can express them 



on \z\ 



1 as 



.4 
B 
B 
A 



(3.8a) 
(3.8b) 
(3.8c) 
(3.8d) 
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We can rewrite the eigenvalue problem ( 13.1 j) in the following equivalent 
forms 



Z- n ^l,n 



z 2 I z 2 Q r . 
R n I 



I Qn 

z~ 2 R n z~ 2 I 



z-( n+1 ^ hn+l 
z-( n+1 ^ 2<n+1 

Z n+1 *>l, n+ l 
Z n+l ^2,n+l 



(3.9a) 
(3.9b) 



Thus, in view of the boundary conditions (13.51) . z n (p n , z n 4> n , z n ip n and z n ip n 
depend on z only through z 2 . Similarly, ( 13. 3 j) can be rewritten as 

[z n+l ^ n+1 ^ +1 $ 2 ,„ +1 ] = [z n $i, n z^ 2 , n ] I" : ~ l :2Q " 



Rn I 



(3.10a) 



[z-(" +1 )$ lin+1 z-( n+1 ^ 2tTl+1 ] = [z- n ^ n z~ n ^ n ] 



I Qn 

z- 2 R n z~ 2 I 



(3.10b) 



so z~~ n ipn d an d z n ip^ d depend on z only through z 2 . Therefore, relations 
(13. 8p imply that the scattering data A, B, B and A are even functions of z 
(cf. [37]); they can be denoted as A(p), B(fj,), B(/i) and A(p), where 

/i = z 2 . 

We introduce the following representations of the Jost solutions ip n and 



rv 

oc 



= n 

rx 

oo 

'<", If 



'ill 






' ' 




' ' 


Ri 


I 




I 




I 



^^i k+l K{n,n + k), (3.11a) 



fc=0 







' I ' 




' I ' 















+ ^ ^T k ~ l K{n 1 n + k), 

(3.11b) 



fc=0 



which are assumed to be uniformly convergent in \fi\ < 1 and > 1, re- 
spectively (cf. (13. 2p ). Here and hereafter, the order of the matrix product is 
defined as 



rv 

m 



m 



Y\ Xi :— X n X n+ i ■ ■ ■ X m , \ J Xi :— X m X m _i ■ ■ ■ X n , m > n, 



7 In the area of orthogonal polynomials, the Ablowitz-Ladik eigenvalue problem in a 
similar rewritten form was studied by G. Baxter in the early 1960s after the seminal work 
of G. Szego, see [HI13]. 
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and the "integral kernels" K(n, m) and K(n, m) are \x- independent 21 x I 
matrices denoted in terms of I x I matrices as 



K(n, m) 



Ki{n, m) 
K 2 (n,m) 



K(n, m) 



Ki{n, m) 
K 2 (n,m) 



m> n. 



We substitute (13.1 lap and ( 13.11b|) into (13.9aj) and (13.9b|) . respectively. Noting 
that they are identities in /i, we can express the "integral kernels" recursively 
in terms of the potentials Q n and R n ; the most important relations are 



Ki(n,n) = Q n , 

oo 

K 2 (n,n) = ^2,RjQ j+ i, 



(3.12) 



j=n 



and 



K 2 {n,n) = R n , 

oo 

Ki(n,n) = y^QjRj+i. 



(3.13) 



]=n 



In a way similar to (13. lip , we can also express the other (adjoint) Jost 
solutions as power series in either /i or fi' 1 . Indeed, noting the identity, 



Rn I 



we obtain 



n-l 



z n </>n= n 

i=— oo 
n-l 

z- n K= n 



—fi 1 R n {I — Q n Rn 



(I- QiRi)- 1 _ 
-Ri (I — QiRi) 



Qn(I RnQn) 




' I 0' 


1 (I -RnQn)' 1 




I 



-\iQi (I - RiQi) 1 
^{I-RiQiY 1 





' I ' 








, (3.14a) 



i=— oo ■- 



fi- 1 (I - QiRi)- 1 -Qiil-R&iY 1 
-fi^Riil-QiRi)' 1 (I -RiQi)' 1 










-I 



(3.14b) 



z - n r n d = [i o]H 

oo 



(I - QiRi)" 
—Ri (I — QiRi 



-\iQi (I - RiQi 



u- 1 (I - QAy 1 -Qi (I - RiQi 



-fi Ri (I — QiRi 



(I — RiQi 



etc. 
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Thus, ( I3.8a|) and (13.8d|) imply that A(fi) and A(fi) can be written explicitly 

as 



oo 



A(n)=[l 0} H 



l=— oo ■- 



oo 



a^)=[o i] n 



l=— oo ■- 



(/ - Q l R i y 1 -/iQi {I - RiQi)- 1 
-Ri (I - QiRiY 1 ^{I-RiQiY 1 



fi- 1 (/ - Q^r 1 -Qi (I - RiQi)' 1 
-^- l R l (I-Q t R t y 1 {I-R l Q i y 1 













(3.15a) 





" ' 




I 



(3.15b) 



Therefore, as long as Q n and R n decay sufficiently rapidly as n — > ±oo, 
z n <fi n and z~ n tfj^ d are analytic on and inside the unit circle < 1), and 
z~ n (p n and z n ipn d are analytic on and outside the unit circle (|/i| > 1). Conse- 
quently, A(a) and A(n) can be analytically continued for \n\ < 1 and > 1, 
respectivelyjj A more precise discussion on the analytical properties of the 
Jost solutions can be made using a discrete analog of the approach in jl]; 
that is, we can rewrite the eigenvalue problem in the form of linear summa- 
tion equations and discuss the convergence of their Liouville-Neumann-type 
series solutions. However, we omit such a discussion in this paper. 



3.3 Gel'fand— Levitan— Marchenko equations 



We multiply (I3.7ap and f!3.7bj) from the right by z n A(fx) 1 and z n A(/i) 1 , 
respectively, to obtain 

[z^MAifi)- 1 = [z n $ n ](fi) + [*-^ n ](/i)S(/i)A(/i)-y\ (3.16a) 
[z~ n (f) n )(iJ l )A(iJ l )~ 1 = -[z-^if,) + [z n y n ](fx)B^)A^)- l ^ n . (3.16b) 

Here, "(/i)" emphasizes that the argument of the functions is fi(= z 2 ) rather 
than z. 

Then, we substitute the summation representations (13. lip into the right- 
hand side of (I3.16al) and operate with 



2m 



(m > n) 



(3.17) 



c 



8 Here, we use the term "analytic continuation" loosely. See appendix B of [44 for a 
rigorous treatment of "analytic continuation" in the delicate case, e.g., when Q n and R n 
do not decay exponentially fast as n — > ±oo. 
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on both sides. Here, C denotes the counterclockwise contour along the unit 
circle \fx\ = 1. Thus, we obtain 

J(n, m 

where 



K(n, m) + 



O 

F c (m) 



+ ^2K{n,n + k)F c {m + k + l), (3.18) 



k=0 



J(n, m) 



< c (mJ 



2vri£ 



^> n ](/i)A( M )-V m - n d/i 



(3.19) 



1 

2vri J c 



<b 5(/i)A(/i)~V m d/i. 



Because of the analyticity of [z n (f) n ](fx) and A{p) in < 1, we can evaluate 
J(n, m) using the residue theorem. Recall that the inverse of the matrix 
A(/x) is given by 

where the tilde denotes the adjugate (i.e., transposed cofactor) matrix. Thus, 
the singularities of the integrand in (I3.19P are determined by the zeros of 
det A(fi). For simplicity, we assume that the matrix function ^(/i) -1 only has 
isolated simple poles in < 1, denoted as {//i,// 2 , • • • , A*at}, and is regular 
on \n\ = ljj In fact, the more general case where A(/z) -1 also has higher 
order poles can be recovered by taking a suitable coalescence limit of two or 
more simple poles afterward. 

In the neighborhood of fx = fij, we can expand A(fx) and ^4(/x) _1 as (cf. jl5l 

) 

A(ji) = Aim) + m)A'{fij) + 0((/i - fij) 2 ), det A^) 
1 



A(v)- 1 



-A 



(-i) 



+ Af + 0(/i - n s ) 



where 



fl — flj 



A0mj)A^ = O, A(fj, j )Af > + A'^Afv = I. 



(o) 



A^ 



0. 

(3.20) 



Thus, using f!3.5bl) . (I3.6bl) and f)3.8ap . we obtain 



n„/,ad 



z n <p n A) 



(-i) 



A(fj)A 



(-i) 



O 



[0 0] at fx 



fX j. 



9 We should not assume such a strong condition as det A(fi) has only simple zeros, which 
was assumed in our previous papers [331123] • Indeed, a zero of multiplicity k of detA(/x) 
may be cancelled by a zero of multiplicity k — 1 of A(/x) to give a simple pole of A(fi)" 1 . 
However, this correction does not affect the validity of the formulas in 23,23]. 



17 



Because z~ n ip^ d consisting of Z rows satisfies the boundary condition in f!3.6bj) 
and the eigenvalue problem f!3.10bj) . the rank of [^ _ri V'n d ](A i i) is equal to / for 
all n G Z. Similarly, the rank of [z n i/j n ](fij) is / for all fi6Z. Therefore, there 



exists an Z x Z matrix Cj such that 



(3.21; 



Here, Cj must be n-independent, because both z n (p n and z n ip n satisfy the 
same eigenvalue problem f!3.9bl) . The matrix Cj can be intuitively considered 
as B(/j,j) lim^^. (/i — fij)A(fi)~ 1 , but it is, in general, different from the naive 
residue lim^^ (// — /i,)£>(/i)v4(/i) _1 . Indeed, can have a discontinuity 

at /i = fij. Because (j) n A ^ ^ and VnCj vanish exponentially for n — >■ — oo 
and n — > +00 respectively, each nonzero column vector of the 21 x Z matrix 
(j) n A^- = tfj n Cj at z 2 = Hj gives a bound state in the potentials Q n and R n . 

Therefore, using the residue theorem with the aid of (I3.20p and (13.211) . 
we can compute the right-hand side of (I3.19P as 



J(n, m) 




O 
F D (m) 



+ t*) +l K{n, n + k)j C jf i? 

K(n, n + k)F D (m + k + 1) 



00 

E 

k=0 



where 



N 

F D (m) : ^O/'-'- 

Substituting this expression for J(n,m) into (13.181) . we obtain a linear sum- 
mation equation of the Gel'fand-Levitan-Marchenko type, 



K(n, m) + 



O 
F(m) 



+ ^ K ( n > n + k )F(m + k + l) 



k=0 



O 

o 



m > n. 
(3.22) 



Here, F(m) is defined as 

F(m) := F c (m) + F D (m) 



1 r - 
— f 5(/i)A(/i)-V m d/i-^a 

J C A 1 



(3.23) 
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Note that F c and F D correspond to the contributions of the continuous and 
discrete spectra, respectively!^! 



Remark. Using the expressions fl3.14a|) and ( I3.15al) . we can evaluate [z n (p n ] (/i) 
and A(fj) in the limit fj, -> as (cf. [1511251130] ) 



lim[z n n ](/i) 



I 

— Rn-l 



n-1 



} \ (i-Q l R t y 1 



oo 



hmAQj) = j f (I - QiRiY 1 . 

i=— oo 

Because lim^o A(fj,) is invertible (cf. (13.41) ). we have 7^ 0, j = 1, 2, . . . , N. 
Note that instead of (13.171) . we can operate with ^ § c d/i on (j3.16aj) with 
(13. lip . Thus, we can express lim M _ ilO [2; n 0„](/i)A(/i) _1 as 



lim^K/i)^^)- 1 



F(n- 1) 



fc=0 



Ki{n, n + k) 
K 2 (n,n + fc) 



F(n + k). 



From the above three relations, we obtain 



I 

—Rn-l 



rx 

00 



H (I - QiRi 



F(n- 1) 



fe=0 



Ki(n, n + k) 
K 2 (n,n + k) 



F(n + k). 



The nonlocal quantities on the left-hand side can be used to transform the 
matrix Ablowitz-Ladik lattice ( 12. 3 p to other systems [7]. 

Next, we substitute the summation representations (13. lip into the right- 
hand side of (13.16bl) and operate with 



27ri£ 



d/i n 



n—m—2 



(m > n) 



(3.24) 



c 



on both sides. Thus, we obtain 



J(n,m) = —K(n,m) + 



F c (m) 
O 



+ ^2K(n,n + k)F c {m + k + l), (3.25) 



k=0 



10 



In this section, we often follow the notation of Ablowitz et al. [6, lMTB]. 
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where 



J(n, m) 



1 

2tti 



[z->„](/i)A(/i)-V n - m - 2 d/. 



(3.26) 



F c m := — 



1 f =, 



— f B^)A^)-^- m - 2 d^. 
J c 



Because of the analyticity of [z n n ](/i) and A(fjL) in \fi\ > 1, we can evaluate 
J(n, m) using the residue theorem. We assume that A(fi)^ 1 only has isolated 
simple poles in \fi\ > 1, denoted as {pi, /Z 2 , • • • , A%h and is regular on = 1. 
In the neighborhood of \i = Jlj, we expand A(p) and A(fi)^ 1 as 



= A(nj) + (fj, - n^A'ifij) + 0((/i - ^) 2 ), det A(^) = 0, 



A(/x) 



-i 



1 



r(o) 



(-i) 



(3.27) 



where 



Thus, using (13.5bj) . (13.6b[) and ( 13.8dj) . we obtain 



n„/,ad 



z n il). 



z~ n <\> n A { r x) z- n i; ri 



-a{h)a\-v o 



= [ O O ] at /i = /x j . 
In the same way as the derivation of (13.2ip . there exists an n-independent 



I x I matrix Cj such that 



(3.28) 



Therefore, using the residue theorem with the aid of (I3.27P and (JH3 
we can compute the right-hand side of (I3.26P as 



N 



■vj-m-2 



J(n,m) = -Y^i^nK^Cjfij 

oo 

k=o J 

oo 

^2 K( n ' n + k)F D (m + k + l) 



3=1 
N 

£ 





" I " 


{ 






— -m-2 



F D (m) 
O 



fc=0 
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where 



N 



F D (m) := y^Cj/i 



— -rn-2 



3=1 



Substituting this expression for J(n, m) into (I3.25p . we obtain another linear 
summation equation of the Gel'fand-Levitan-Marchenko type, 



Kin, m) 



F(m) 
O 



+ ^ K{n, n + k)F(m + k + l) 



k=0 



O 

o 



m > n. 
(3.29) 



Here, F(m) is defined as 

F(m) := F c (m) + F D (m) 



It-- N 

— f B^)A^)-^- m - 2 d^ + Y J C 3 Jl, 

71 ^ C 1 = 1 



-m-2 



(3.30) 



Remark. Using the expressions (13.14bj) and (13.15bj) . we can evaluate [z ™0„](/i) 
and A(fi) in the limit \fi\ — > oo as (cf. [40J) 



lim [z n (j) n ](n) 

»oo 



Qn-1 
-I 



n-1 



H {I-FUQiY 1 



oo 



lim A(fj) = FT (I-RiQi)' 1 . 

i=—oo 

Because lim^i^oo A(fj,) is invertible (cf. (I3.4p ). A(fi)^ 1 does not have poles at 

infinity, which is indeed consistent with the computation for J(n,m). Note 
that instead of (13.241) . we can operate with ^r^d/i/x -1 on (13.16bj) with 

(13. lip . Thus, we can express \xm\ tl \^. 00 [z~ n <p n )(iJ l )A(iJ l )'' 1 as 



lim [z-^MA^Y 

\ft\— >oo 



-F(n- 1) 
I 



From the above three relations, we obtain 

rx 

oo 

n ( j - RiQt) = 



—Qn-l 
I 



-F(n-1) 
I 



fc=0 



k=0 



K\(n, n + k) 
K 2 (n,n + k) 



Ki{n, n + k) 
K 2 (n,n + k) 



F(n + k). 



F(n + k). 
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The nonlocal quantities on the left-hand side can be used to transform the 
matrix Ablowitz-Ladik lattice ( 12. 3 p to other systems [7]. 



The Gel'fand-Levitan-Marchenko equations (I3.22p and (13.291) relate the 
scattering data to the potentials Q n and R n through (I3.12p and (I3.13p . The 
required set of the scattering data is given by -B(/i)A(/i)~ 1 and £>(/i)A(/i) _1 
for = 1, {fij, Cj}j=i t 2,...,N and {fij,Cj}j =12 wn i cn define F{m) and 

F(m) as in (13.231) and (13.301) . Because F(m) and F(m) are "linear" in the 
scattering data, we can consider a linear superposition of different sets of 
scattering data at this level. In fact, we will show in the next subsection that 
F(m) and F(m) satisfy linear evolution equations. 



3.4 Time evolution 

Under the rapidly decaying boundary conditions ( 13. 2p . the temporal part 
of the Lax representation (12 .2p for the Ablowitz-Ladik lattice (12. 3p has the 
asymptotic behavior: 



" *l,n " 




. ^2,n _ 


t 



-li + l)bl O 
O {l-n- x )al 



*2,n 



as n — > ±oo. 



This can be used to fix the time dependence of the leading order terms in \l/i jn 
and ^2,n (cf. (13. 5p ). Thus, we can introduce the explicitly time- dependent 
Jost solutions cpn and (f>$ as 



' I 

O 



and and ipn^ as 

z n ^ := e^ +1 ) u z n ^ n -»■ e ( -^ +1)6 * 



(l-/4 _1 )at 



o 
-I 



o 
I 



> as n — > — oo (3.31a) 



7 
O 



> as n — > H-oo, (3.31b) 



respectively; they satisfy both the eigenvalue problem ( 12. ip and the time- 
evolution equation 
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To determine the time dependence of the scattering data, we rewrite the 
defining relations ( 13. 7p as 

Note that f% <pn\ ip® and satisfy the same equation (12. 2p and the 
columns of ij)^ and tpn are linearly independent. Thus, the time depen- 
dences of A, B and A, B for = 1 are given by 

A(ji, t) = A(ji, 0), B(ji, t) = B(ji, ) e [(M-i)''+(i-M- 1 Ht (3.32) 

and 

A(jjl, t) = A(p, 0), B{fi, t) = B(ji, ) e -[^- 1 ) 6 +( 1 -^" 1 ) a ]*, (3.33) 

respectively. This implies that A(ji,t), A(ji,t), B(ji,t)e-lb t - 1 M i --i*- 1 )*] t and 
B(fi, t)e^ _1 ) 6+ ^ 1_// ^ are generating functions of the integrals of motion. 

Because the "analytic continuation" of A([a) and A(fj) into the regions 
|/i| < 1 and > 1, respectively, is unique and remains time- independent 
(cf. (I3.15P ). the positions of the simple poles of A^fi) -1 and A(fi)" 1 and the 
corresponding residues, Aj }j =12 N and {/Zj, ^}-_ 12 ^, are also 
time-independent. For (13.211) and (13.281) . we can apply a similar discussion 
as used to obtain (I3.32p and (I3.33p . so the time dependences of Cj and Cj 
are given by 

Cj{t) = C i (0)e^- 1)6+(1 -^ 7l ) al * (3.34) 

and 

Cj(t) = ^(0)e- [( ^- 1)6+(1 -^ 1)a] *, (3.35) 

respectively. 

Substituting (l3T32|) - fl3T35l ) into ( 13T23|) and ( 13T30|) . we obtain the explicitly 
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time-dependent forms of F{n) and F{n) as 

F(n,t) = — l S( A t,0)A(//,0)-V n e [(At " 1)6+(1 "^ 1)ol *d^ 
2vri Jc 

N 

~ Yl Q(0)/i"e [( ^" 1),,+(1 "^ 1)a] *, (3.36) 

3=1 

F(n,t) = — / 5(/i,0)l(/i,0)-V" n " 2 e" [( ^ 1){,+(1 ^ 1)al *d/i 
2vri J c 

TV 

+ S ^(0)/7^-V^- 1)6+(1 ^ rl)a] *. (3.37) 

3=1 

Thus, it is easy to see that F(n, t) and F(n, t) satisfy the pair of uncoupled 
linear evolution equations: 

dF(n t) 

1 ' ; - bF(n + 1, t) + aF(n -l,t) + (b- a)F{n, t) = O, (3.38a) 



dt 
dF(n, t) 
dt 



aF{n + 1, t) + bF{n - 1, t) + (a - 6)F(n, t) = O. (3.38b) 



Note that these equations coincide with the linear part of the equations for 
R n and Q n (see ( 12. 3p ). In addition, F(n,t) and F(n,t) are required to de- 
cay rapidly asn-y +oo so that the Gel'fand-Levitan-Marchenko equations 
(13^21 and (EEZg) are well-posed. 

Because of the linear nature of the sum terms in (I3.36P and (13.371) . we can 
take a coalescence limit of two or more simple poles of A(n)~ l and A(fi)~ l 
directly. Thus, we obtain the following generalized expressions for F(n, t) 
and F(n,t): 

F(n,t) = —i B(^,0)A(fi,0)~ 1 fi n e [ ^ 1)b+{1 ~^ 1)a]t dfi 
2tti Jo 

N Mj / a \ k 

n [(^-l)b+(l-^" 1 H* 



iv / a \ K 

EE c f (") (£} & 

j=i fc=o \ r-3 / 



F(n,t) = ^-j> 0)^^, 0)-V~ n-2 e -[( ^ 1)6+(1_M 1)a] *d// 



>c 

N L 



+ EE^c)(4)f. 

j = l fc=0 \ f"3 / 



3 ^ ' 
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These expressions encompass the most general case where and A(fi)^ 1 

have arbitrarily higher order poles. Moreover, they satisfy the same linear 
equations (I3.38|) as the original expressions (I3.36P and (I3.37p . 

Instead of using the partial differentiation with respect to fij and Jlj, one 
can consider the matrix functions x n e bt ( x ~ I * )+at( - I ~ x > and Y~ n ~ 2 e~ bt ( Y ~ I >~ at ( I ~ Y ' 
with constant invertible matrices X and Y in Jordan normal form. Indeed, 
they satisfy the linear equations of the form (I3.38P , so linear combinations of 
the independent elements of each matrix function can be used to replace the 
sum terms in F(n,t) and F(n,t). Readers interested in such an approach 
are referred, e.g., to |^f ]H71 - l52] . 

3.5 Exact linearization 

To reconstruct the potentials Q n and R n from the scattering data through 
(13.1 2p and (I3.13p . we rewrite the Gel'fand-Levitan-Marchenko equations 
(I3.22p and (I3.29P as "closed" linear summation equations for Ki(n,m) and 
K 2 (n, m). Thus, the general solution formulas for the matrix Ablowitz-Ladik 
lattice f)2.3p can be presented in the form: 

Q n = K 1 (n,n), (3.39a) 
R n = K 2 {n,n), (3.39b) 

oo oo 

Ki(n, m) = F{m) — K\{n, n + i)F{n + i + k + l)F{m + k + 1), m > n, 

i=0 k=0 

(3.39c) 

oo oo 

Kiin, m) = —F{m) — K 2 {n, n + i)F{n + i + k + l)F(m + k + 1), m > n. 

t=0 k=0 

(3.39d) 

Here, the time dependence of the functions is suppressed; F{n) and F{n) 
are solutions of the linear uncoupled system (13.381) and decay rapidly as 
n — > +oo. More generally, the set of formulas (I3.39P can provide the solutions 
for any flow of the matrix Ablowitz-Ladik hierarchy if F(n) and F(n) satisfy 
the linear part of the equations for R n and Q n , instead of (I3.38p . Hence, 
(I3.39P realizes an exact linearization of the matrix Ablowitz-Ladik hierarchy 
in the sense of [TD] (also see [SUED and Proposition IA.2[) . As long as F(n) 
and F(n) decay rapidly as n — > +oo, so do Q n and R n determined by (I3.39p . 
However, the requirement that Q n and R n should also decay as n — > — oo 
imposes nontrivial conditions on F(n) and F(n), which will be touched upon 
in the next subsection. 
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3.6 Multisoliton solutions 



To construct exact solutions in explicit form, we consider the special case 
of B(fi) = B(fi) = O on = 1; this is preserved under the time evolu- 
tion (cf. (I3.32p and (13.331) ) and corresponds to the reflectionless potentials 
(cf. (13. 7p ). Moreover, we assume that A{h) 1 and A^) 1 only have simple 
poles (see [4T] for the more general case). Thus, we can set 



N N 

F(n,t) = F(n,t) = ^C,(t)/Z7- 

3=1 3=1 



(3.40) 



where the time dependences of Cj and Cj are given by (I3.34p and (I3.35P . We 
also set 



iV 



N 



K x {n,m\t) = J2 G j( n ^)N m ~^ K 2 (n,m;t) = ^tf^n,*)//™ (3.41) 

3=1 3=1 

and substitute all these expressions into (I3.39cp and ( 13.39d[) ; recalling that 
\fij\ < 1 (j — 1, 2, . . . , N) and \JIj\ > 1 (j = 1, 2, . . . , N), we can evaluate the 
infinite sum. Thus, we obtain a linear algebraic system for determining Gj 
and that for Hj as 



Gi/i 1 n 2 G2H2 n 2 ' ' 

r~< —-n~2 ri — — n-2 



IN 



u 



Nl 



u 



NN 



N^ N 



—-n-2 



(3.42a) 



and 



V n ■■■ V 1N 
Vjvi • • • Vnn 



[ Cxji n C2H2 "' Cn^n ] ■ 
(3.42b) 



Here, all the entries in (I3.42p are I x I matrices; the block matrices U = {Uji 
and V = {yjk)\<j,k<N are defined as 



l<j,k<N 



N 



i=i ( 1 _ ^ ) ( 1 _ £i 



to 



■Ci{t)C k (t) 
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and 



N 



V jk :— 5 jk I - 2j 



i=i 1 _ 



1 - 



Ci(t)C h (t), 



respectively. Here, Sjk denotes the Kronecker delta. Thus, using f)3.39al) . 
fl3~59b|) and (jS3H) , we obtain 

Q„(t) = ifi(n,n;t) 



I 
I 



Cx{t)K 



-n-2 



u 



11 



%1 



IN 



-1 r 



(3.43a) 



R n {t) = K 2 {n,n; t) 

— [Hi^i ■ ■ ■ H N fi N ] 



[C^rt ■■■ C N (t)fJL%] 



11 



IN 



Vm ■■■ V T 



NN 



-1 r 



(3.43b) 



This provides the multisoliton solutions of the nonreduced matrix Ablowitz- 

Ladik lattice (I2.3P ; some additional conditions on {fij, Cj}j=i i2 ,...,N and {Jlj, Cj} . =1 2 ^ 

need to be satisfied for (I3.43P to exhibit solitonic behavior. 

In the simplest nontrivial case of N = N = 1, we obtain the one-soliton 
solution of (12.31) in the form: 



Q n (t) = D{n,t) 



R n (t) = D(n,t) 



-1 -1 



I- 



mi 



1 _ ffi. 



Hi 



; D(n,t)D(n,t) 



-1 -1 



1 _ 



; D(n,t)D(n,t) 



(3.44a) 



(3.44b) 
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where D(n,t) := dW/Zr^V^-W 1 -^ X W and D(n, t) := C 1 (0)^e^-^ b+( - 1 -' l i 1 ^ t . 
For this solution to decay also as n — > — oo, we require that lim^.oo Q n (t)l = 
for any n-independent column vector / of dimension I and similar for R n (t). 
Thus, considering the Maclaurin series for (/ — X) -1 where X is the corre- 
sponding matrix in (I3.44a[) or (13.44b|) . we obtain the following conditions on 
the kernels of the / x I matrices C\ and C\. 



and 



Ker (ddd) = Ker (d) 



Ker ddd = Ker(d 



Note that these conditions remain invariant under the time evolution (cf. ( I3.34p 
and (13.351) ). They can also be written in a more easy-to-understand form: 



Ker (Ci) n Im d = Ker Ci n Im (d) = {0} 



Consequently, we have rank ( d ) = rank (C 



For general values of N and N, it is rather difficult to grasp the condition 
that Q n and R n given by (I3.43P should also decay as n — > — oo. Thus, we 
take a different route. In view of the first component of (I3.29P and the second 
component of (13.221) . relations (13.401) and (I3.4ip together with (13. lip imply 
that 



— -ra-2 



Gj(n,t)jl- n - 2 



1,2,. ..,N, 



1,2,. ..,N. 



Thus, Gj and Hj are closely related to the bound-state eigenf unctions. Ow- 
ing to the connection formulas (I3.28P and ( I3.2ip as well as the boundary 
conditions ( 13. 5 aft . GjJiJ n ~ 2 and Hjfx™ must decay as n — > — oo. If this is sat- 
isfied, then Q n and R n in ( I3.43P naturally vanish asn-y — oo. The relations 
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(I3.42p for determining Gjfij n and Hj/j,™ can be rewritten as 



[G 1 G 2 ■ ■ ■ ] - [ Gi G 2 ■■ ■ G^ ] 



n+l n-2 

Ml Ml T 

I n 1 

mi 



n + l n-2 

Mi u_ 
1 AT 7" 

Ml ^ 



n+l n-2 -i 

Mat Mi 7- n 

Ml 



n + l n-2 

^ ^ r 



x 



and 



mi— mi - 1 - 1 



M^-Mi AT 



1 /^r /~i 1 /m /m 
=— O/vrOi • • • — <~ / A r< ~ / 



HI— UN 



C*i C*2 • • • 



(3.45a) 



[ i^i i7 2 • • ■ iljv ] — [ -Hi #2 ' " " -Hjv ] 



nTT-n-3 

Ml Mi j 



Ml M- 



1-4 



''i 



MjyMi r 
i _ iiiv J 



M — 

N 



,,n TJ-™- 3 



X 



Ml Mi 



i 1 _i C\Cn 

Mjv Mi 

i 1 _i C^Cn 



[ d C 2 ■ ■ ■ C N ] . (3.45b) 



Then, we multiply both sides of fl3.45a[) from the right by an n-independent 
column vector of dimension I x N and consider the limit n — > — oo. Thus, we 
obtain the condition 



Ker 



mi— mi 



dd 



-C X C; 



M^T-Mi 1 N 



1 /~i /~t 1 /~i /~t 

Mi-MJV JV 1 M^-Miv IV JV 



C Ker 



d C 2 ■■■ C 



N 



Similarly, from (13.45bj) . we obtain 



Ker 



Ml Ml 



Ml M- 

JV 



i 1 _i CiCn 

MJV Ml 



1 1 _1 CjyCjV 

fjv M — 

JV 



(3.46a) 

C Ker [ d C 2 • • • C N ] . 

(3.46b) 
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The two conditions (I3.46ap and (13.46b[) can be combined and simplified to 
provide more easy-to-understand conditions on the scattering data in the 
reflectionless case. For this purpose, we need the following lemma. 

Lemma 3.1. Define the (M — 1) x M matrix elements dy G C as 
(cii, bj) 



:- . 

A; — V a 



i G 0*i, i 2 , • • .,«m-i}, j G {ji,h, ■ ■ -,3m}- 



Here, Aj and Uj are parameters, and bj are nonzero vectors of dimension 
I and ( • , • ) stands for the scalar product. For d^ to be well-defined, we as- 
sume Aj 7^ Vj for all i and j, but we do not require X ia ^ X ij3 or Uj a ^ Vj„ 
for a 7^ p. Instead, we assume the following condition: for any subset 
{h, k 2 ,..., k^} C {ji, j 2 , • • -,3m} such that u kl = u k2 = ■ ■ ■ = v^, the vec- 
tors 

bk 1} b k2} . . . , b ki 

are linearly independent. Then, if the equality 



M 



Q=l 



d 



*U1 



dilja-l di 1 j a + 1 



d 



diM-iji ' ' ' di M _ 1 j a _ 1 di M _ 1 j a+1 ■ ■ ■ di M _ x j M 



K = o 



is valid, all the scalar coefficients must be zero, where | ■ | stands for the 
determinant. In other words, the above vector equation holds true only in 
the trivial case; note that this equation can be written compactly as 



di 
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h 3M 

diijM 



d; 



di 



l " l M-l3l l M-ljM 

using the Laplace expansion formally. 



We omit the proof of this lemma. To obtain useful information from 
the conditions (I3.46p . we first remove the trivial subspace of the kernels 
commonly contained on both sides. From a given I x I matrix W, we extract 
the maximum number of linearly independent column vectors to form an 
I x rank(V*/) matrix W^ c \ Similarly, we extract the maximum number of 
linearly independent row vectors from W to form a rank(W) x I matrix . 
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With this notation, (I3.46a|) and f!3.46bj) can be rewritten in more compact 
forms as 



Ker 



and 



Ker 



" 1 c^C^ 
mi -mi 1 1 

i c^C^ 

Ml -MiV N 1 
_L L I °1 

n mi 

— -J— jy 1 



1 TV 



C C 





'Ml 


1 






MiV 


1 




1 


j_ 


MAT 


Ml 


1 




_J 


l 


MJV 


~ — 

N 



N 



C Ker 



C 



(_c) 

iV 



(3.47a) 



'1 °7V 



''A' 



C Ker 



c 



(c) 
JV 



(3.47b) 

With the aid of Lemma 13.14 we can prove the following two propositions. 
Proposition 3.2. Assume that 

N N 

rank (Cj) < rank ^Cjj > 



3=1 



and the condition (I3.47a|) is satisfied. Then, the above inequality becomes 
an equality, 

N N 

rank (Cj) = rank (Q) > ( 3 - 48 ) 



3=1 3=1 

and the matrix on the left-hand side of f!3.47ap must be invertible, i.e. 



i c^C^ ■ ■ ■ i c^C^~ 

Mi-Mi 1 1 M j\7 Ml 1 N 



1 c^C^ 

Ml -Mat ^ 1 



M^-MJV 



N 



^0. 



(3.49) 



Proposition 3.3. Assume that 



N N 

rank (Cj) > rank (pj) > 

3=1 3=1 
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and the condition f!3.47bj) is satisfied. Then, the above inequality becomes 
an equality, 

N N 

^rank(0,) = J^rank (pjj , 

3=1 3=1 

and the matrix on the left-hand side of (13.47bj) must be invertible, i.e. 



l r<(*)n{ c ) 
_i r ^i ^i 

Ml Mi 



— — Jv 1 

Ml M- 

AT 



i n( r )r<(. c ) 

1 1 ^1 



MJV Ml 



i r'Wr'W 

^^ U N N 



(3.50) 



By combining Propositions 13.21 and 13.31 we arrive at the following theo- 
rem. 

Theorem 3.4. Assume that (I3.43P provides the multisoliton solutions of the 
matrix Ablowitz-Ladik lattice ( 12. 3p . which decay as n — > ±oo and produce 
the bound states of the associated eigenvalue problem (13. ip . Then, the scat- 
tering data must satisfy the three conditions ( I3.48p -( 13.5(jp . 

Note that the time evolution does not change these conditions (cf. ( I3.34p and 
(J335J)). 

Remark. Here, we only considered the case where the time variable t is 
fixed at some finite value. In the limits t — > ±oo, the solitons generally sep- 
arate from one another and restore their original shapes. Thus, for ( I3.43P 
to describe proper multisoliton collisions through the passage of time, we 
have to impose additional conditions, i.e., the above conditions for subsets of 
{fij, Cj}j = i t 2,... t N an d {yUj, Cj}. =l 2 jy. Some relevant results were obtained 
independently in [4"Tj . 

3.7 Complex conjugation reduction 

When b = a*, the matrix Ablowitz-Ladik lattice (I2.3P allows the complex 
conjugation reduction R n = aQ* n with a real constant a (cf. |7J). In particu- 
lar, the simplest reduction R n = —Q* n can be realized in formulas (I3.39P by 
identifying F(n) with the complex conjugate of F(n), i.e. 

F(n) = {F(n)Y, (3.51) 
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which is naturally preserved under the time evolution ( I3.38P with b = a*. 
Indeed, this relation can be derived by exploiting the symmetry of the eigen- 
value problem (13. ip with R n = —Q* n ] that is, if 

" *l,n(z) " 

is an eigenfunction, then 

, r -*2,»(i/**) r 

L *l,n(lA*) 

gives another eigenfunction of the same problem. Thus, we can reflect this 
symmetry in the Jost solutions and the scattering data to confirm (I3.5ip . 

In particular, the iV-soliton solution of the matrix Ablowitz-Ladik equa- 
tion, 

Qn,t - aQn+i + a*Qn-\ + (a - a*)Q n - aQ n Q* n Q n+1 + a*Q n ^Q* n Q n = O, 
is obtained by setting b = a*, N = N and 

V>j = A. CjMW 2 = ~ {Q(*)>* » J = i, 2, . . • , iV 

in formula fl3.43a|) with (I3.34p . The reduced set of scattering data is required 
to satisfy the conditions ( I3.49P and (j3.50P ; in fact, they are equivalent under 
this reduction. In addition, ( I3.49P (or (I3.50p ) for subsets of {fij, Cj}j=i t 2,...,N 
should also be satisfied. Throughout this paper, we do not discuss the issue 
of regularity of solutions and the term "soliton solution" is used in a broad 
sense. That is, it may have singularities at some values of the independent 
variables n and t. 

4 Solution formulas for the derivative NLS 
lattices 

4.1 Solutions of the space-discrete Gerdjikov— Ivanov 
system 

In this subsection, we solve the space-discrete Gerdjikov-Ivanov system de- 
rived in subsection 12.21 by using the results in section [3J Note that the 
nonzero parameter \x in the space-discrete Gerdjikov-Ivanov system (12. 5p is 
nonessential; it can be fixed at any nonzero value, say 1, using a simple point 
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transformation and rescalings of the parameters a and h (cf. [ITIED])- in 
addition, as is clear from the defining relation of the Miura map (j2.4a|) . the 
limit // — > is trivial and need not be considered separately. Thus, in the 
following, we consider the space-discrete Gerdjikov-Ivanov system (12.51) with 
\i = 1: 

Qn,t — aQn+i + bQ n ^i + (a — b)Q n + aQ n (P n — P n +l) Qn+l 

— bQ n -i {P n — P n +l) Qn + aQnPnQnPn+lQn+1 — bQ n -iP n Q n P n+ \Q n = O, (4.1a) 

P n j - bP n+1 + aP„_i + (b - a)P n - bP n (Qn-i - Qn) P n +i 

K + ClPn-l (Qn-1 — Qn) Pn + bP n Q n -iP n Q n P n+ i — aP n -\Qn-\PnQnPn = O. (4.1b) 

In section [3j we developed the inverse scattering method associated with 
the matrix Ablowitz-Ladik eigenvalue problem (13.11) under the vanishing 
boundary conditions on the potentials Q n and R n (cf. (13. 2p ). The remaining 
unknown P n in (14. ip can be determined from a linear eigenfunction through 
the simple formula 

Because there is some arbitrariness in choosing the linear eigenfunction, we 
need to specify boundary conditions to determine P n uniquely. Thus, we 
assume that not only Q n but also P n decays rapidly as n — > ±oo: 

lim Q n = lim P n = O. (4.3) 

n— >±oo n— ^±oo 

This is consistent if we choose the linear eigenfunction appearing in the right- 
hand side of (14. 2 p as 



2.ii 



4>n, 



and assume that the scattering data vanishes at // = 1, i.e., B(l) = O 
(see (I3.5P and (I3.7aj) ). In fact, the Jost solution ip n as well as n does not 
satisfy the time part of the Lax representation (12.21) . so it is more appropri- 
ate to use the explicitly time-dependent Jost solutions introduced in subsec- 
tion [331 However, the overall multiplicative factor e(-i M + 1 ) bt as introduced in 
(I3.3ip plays no role in formula (14. 2 p , so in view of (13.1 lb[) , we can express P n 
as 

{oo "\ f oo 

J2K 2 (n,n + k) \ I I + Y,Ki(n,n + k) 
k=0 ) I k=0 

This expression enables us to determine P n from the set of scattering data 
with the aid of the Gel'fand-Levitan-Marchenko equations (13.221) and (13.291) ; 
however, for later convenience, we take an alternative approach. 
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Because the space-discrete Gerdjikov-Ivanov system (14. Xp is "symmetric" 
with respect to Q n and P n , there must be a formula for expressing Q n in a 
manner similar to ( 14 . 2 p . Such a formula can be established by identifying an 
appropriate Ablowitz-Ladik eigenvalue problem that is gauge equivalent to 
the original problem (13. IP ; the corresponding gauge transformation is often 
referred to as a Backlund-Darboux transformation. Then, in the new gauge, 
the roles of Q n and P n are swapped and P n appears directly as a potential in 
the Ablowitz-Ladik eigenvalue problem. In other words, there exists another 
Miura map from the space-discrete Gerdjikov-Ivanov system (14.11) to the 
Ablowitz-Ladik lattice in the form (Q n ,P n ) H- (Q n ,P n )- 

Let us consider the original Ablowitz-Ladik eigenvalue problem (13. ip with 
R n = P n - P n+1 + P n Q n P n+l (cf. fl£HD): 



(AL1) : 



Zl zQ n 
Z 1 (P n — Pn+1 + PnQnPn+l) Z l I 



(4.4) 



Indeed, this can be rewritten as another Ablowitz-Ladik eigenvalue problem: 
(AL2) : 







zl 


Z { — Qn + Qn+1 + QnPn+lQn+l) 






$2,n _ 




Z l P n +l 


z-H 




$2,n+l 



using the gauge transformation defined as 



<3>1 









1) * 1,„ -Qn{I~ PnQn)~ l (*2,„ " Z~ 2 PA,r. 
(I - PnQnY 1 (^2,n- Z~ 2 P n Vi, n ) 



(4.5) 



(4.6) 



9r, 



The explicit form of the transformation matrix g n is unimportant; only its 
asymptotic behavior asn-y ±oo is needed: 



lim g n 

n—>±oo 



1)1 o 



o 



For the original Jost solutions for (AL1) defined as (I3.5p . the Jost solutions 
for (AL2) are given as 



0i AL2) (-) 



z- 2 - 1 



9n^ U \z), <P n A ^(z)=g n( p n ^\z) 

1 



z 



'0. 



(AL2) , 



z 



-2 



-9n^, 



(AL1), 



Z 



(4.7a) 
(4.7b) 
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Indeed, they satisfy both the eigenvalue problem ( 14. 5p and the boundary 
conditions ( 13. 5p . The case of z 2 = 1 can be understood in the corresponding 
limit. 

In view of the aforementioned condition B(l) — O, it is natural to modify 
the defining relations (13. 7p of the scattering data for (AL1) on = 1 as 

0i AL1) = tt ALl) A(v) + ^i AL1 V 1 " 1)B(ia), (4.8a) 
gALi) = $jfu)B{p) - ^ AL1) A(/i). (4.8b) 

Then, relations (14 .7p between the Jost solutions for (AL1) and those for 
(AL2) imply that the defining relations of the scattering data for (AL2) on 
= 1 become 

^) = ^AL 2 ) A(/i)+ ^(AL 2)jB(/i)) (4ga) 
^(AL2) = ^(AL2) (// -1 _ 1)5(A|) _ ^W)^). (4.9b) 

The bound-state eigenvalues are determined by the positions of the simple 
poles of A(fi)~ l in |/i| < 1 and A(fi)^ 1 in | | > 1; the more general case of 
higher order poles can be recovered by taking a suitable coalescence limit. 
Because of the uniqueness of the "analytic continuation", the bound-state 
eigenvalues are common to (AL1) and (AL2): 

/4 AL1) = ^f L2) = to J = 1.2,..., N, 

_(AL1) = _(AL 2 ) = _ ^ J = 1)2 _.^. 

Owing to ( 14. 7p . the corresponding matrices Cj and Cj (cf. (I3.2ip and (I3.28P ) 
for (AL1) and (AL2) can be expressed as 

Cf L1) = fc 1 - 1) Cj, C] AL2) =Q, j = l,2,...,A, (4.10a) 

C/ AL1 ) = Q, Cf L2) = (/J/ 1 - 1) Q, j = 1,2,..., iV. (4.10b) 

By combining the above relations, the functions F(m) and F(m) for (AL1) 
and (AL2) (cf. fl3T23|) and (jS20D ) can be written as 

F (AL1) (m) = F(m - 1) - F(m), F (AL2) (m) = F(m), (4.11a) 

F (AL1) (m) = F(m), F (AL2) (m) = F(m + 1) - F(m), (4.11b) 

in terms of the original F(m) and F(m) defined as (I3.23P and (13.301) . Clearly, 
the modification of the scattering data for (AL1) and (AL2) as described 
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above does not change their time dependences given by (I3.32l) - ()3.35p . Thus, 
each of the pairs (F, F), (F^ AL1 \F^ AL1 ^) and (_p( AL2 ) ; _p( AL2 )) satisfies the 
same linear evolutionary system (13.38 p . 

We can construct the Gel'fand-Levitan-Marchenko equations for the space- 
discrete Gerdjikov-Ivanov system (14. ip by combining (I3.39aj) and (I3.39c[) for 
(AL1) and (E09b]) and ( l3~3Mj) for (AL2), i.e. 

Q n = K[ AL1 \n,n), 
P n+1 = Rt L2 \n,n), 

oo oo 

K( AL1) (n, m) = F (AL1) (m)-J2Yl ( n > n + *)^ (AL1) (n + i + k + l) 

i=0 k=0 

x F (AL1) (m + k+ 1), m > n, 



oo oo 



Kt L2) (n,m) = -F^ AL2 \m) ^'V, n + i)F^ Ah2 \n + i + fc + 1) 

t=0 fc=0 

xF (AL2) (m + Hl), m>n. 

Substituting (14. lip and changing the notation slightly, we obtain 

Q n = JtT(n,n), (4.12a) 
P n = ^(n,n), (4.12b) 

oo oo 

Jf{n, m) = F(m) + ^( n ' n + ^ f^ 71 + * + k + ^ ~ F ( n + i + k )) 

1=0 k=0 

xF(m + Hl), m > n, (4.12c) 

oo oo 

X{n, m) = -F(m - 1) - S ^^ n ' n + ^ { F ( n + i + k + I) - F(n + i + k) } 

i=0 fc=0 

xF(m + k), m>n. (4.12d) 

Note that F and F satisfy the linear part of the equations for P n and Q n , 
which is (13.381) for the space-discrete Gerdjikov-Ivanov system (14. ip . 

In the case of B{ji) = B{ji) = O on = 1, which corresponds to the 
reflectionless potentials for both (AL1) and (AL2), we can solve the Gel'fand- 
Levitan-Marchenko equations (I4.12p to obtain the soliton solutions in closed 
form. The derivation is essentially the same as in the Ablowitz-Ladik case 
described in subsection I3.6| naturally, the multisoliton solutions of the space- 
discrete Gerdjikov-Ivanov system (14. ip can be obtained directly by applying 
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the correspondence relations (14.101) to (I3.43p . i.e. 



Q n (t) = C^Wx 



, rr-n-2 



■ ^ 



-i -i 



IN 



// 



Nl 



NN 



I 

(4.13a) 



P n (t)= [Cxit)^' 1 ■■■ C N {t)ii n N - x ] 



y u ... 
Vm ■■■ 



IN 



NN 



(4.13b) 



Here, all the entries in (I4.13P are I x I matrices; the block matrices ty£ = {f^jh) x< ^ 
and "V = {'Vjk)i<j,k<N are defined as 



%k := S jk I - 



v Ii - in) MK n ~ z 



i=i l 



■Ci(t)C k (t), 



N (i _ 77-11 rr n - 2 n n 



i=i i 



and the time dependences of Cj and Cj are given by (I3.34p and (I3.35p . Note 
that the three conditions (I3.48l) - (l3.50p must be satisfied for (14.131) to describe 
proper multisoliton solutions decaying as n — > ±oo (cf. Theorem 13.41) . In 
addition, we need to impose similar conditions for subsets of the soliton 
parameters so that the solitons interact with each other properly throughout 
the time evolution. 

When b = a*, the space-discrete Gerdjikov-Ivanov system (14. ip allows the 
complex conjugation reduction P n = ic<5*_ 1 / 2 with a real constant o |20j. 
That is, two originally uncoupled systems, (14. ip with nfZ and (14. ip with 
hGZ+1/2, can be related by this reduction to give a single equation with 
n G Z/2. Clearly, the value of a is nonessential, so we set a = 1 and consider 
the reduction P n = iQ*_]/ 2 - This reduction can be realized at the level of 
formulas (I4.12p by setting 

F(n) = -i{F(n -§)}*, 
which is consistent with the time evolution (I3.38P with b = a*. In particular, 
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the iV-soliton solution of the space-discrete Gerdjikov-Ivanov equation, 
Q n ,t - aQn+i + a*Qn-i + (a - a*)Q n + iaQ n \Q*_i - Qn+±) Q^+i 

— lCl*Qn-l ~ Qn+^j Qn ~ a QnQ*_iQnQ* + iQn+l + 0-*Qn-lQ*_ l QnQ* + l Qn = O, 



is obtained by setting b = a*, N = N and 

1 ( 



in formula (I4.13ap with (I3.34j) . The imaginary unit (roman i) should not be 
confused with the index of summation (italic i) in the definition of %jk- The 
reduced set of scattering data is required to satisfy the condition (I3.49[) (or 
f)3.50p ) and its smaller versions corresponding to subsets of the solitons. 

4.2 Solutions of the space-discrete Kaup— Newell sys- 
tem 

In this subsection, we solve the space-discrete Kaup-Newell system ( I2.17P 
derived in subsection 12.31 by applying the results in section [3] Because the 
parameter \i is nonessential in (I2.17p . we set fi(= z 2 ) = 1 and consider the 
space-discrete Kaup-Newell system in the form: 

q n>t - A+ [a (I - q n r n y l q n + b(I + g n _ir n ) -1 q n -i] = O, (4.14a) 
r n>t - A+ [b (I + r n g n _i) _1 r n + a(I - 7V-ig„_i) _1 r n _J = O. (4.14b) 

Recall that A+ denotes the forward difference operator: A^/ n := f n+1 — f n . 
We assume vanishing boundary conditions at spatial infinity: 

lim q n = lim r n = O. (4-15) 

n— »±oo n— >±oo 

Propositions 12.11 and 12.21 imply that the solution of (14.1 4p can be obtained as 
a - A+ (V 1)- V 2)> 1 r -f^V'-^-'^V 1 

(/« ~ a n *l,n ^l,n 5 'n — \ ^2,n ^ 2,n ^l,n ^l,n ) 

(4.16) 

using two linearly independent eigenfunctions of the linear problem, (12. ip and 
( 12.21) . associated with the Ablowitz-Ladik lattice (cf. (12.61) ): this is in contrast 
to the space-discrete Gerdjikov-Ivanov system, which can be derived using 
only one linear eigenfunction as described in subsection 12.21 Proposition 12.31 
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implies that r n in (I4.16P can be rewritten in the difference form as q n ; to 
see this explicitly, we need to identify an appropriate linear problem for the 
Ablowitz-Ladik lattice, which is gauge equivalent to the original problem. 
It turns out that the gauge transformation connecting (AL1) and (AL2) 
considered in subsection 14.11 plays the desired role. 

Suppose that the two eigenfunctions appearing in (14.1 6p satisfy (AL1). 
Moreover, we choose the linear eigenfunction appearing in (14.21) as 



p = ^^~ 1 



M (=* 2 )=l 



Recall that (AL1) and (AL2) involve this P n and are connected through the 
gauge transformation (14. 6p . Thus, the two linear eigenfunctions for (AL2) 
can be introduced as 



1,71 



(1) 
2,71 



(1) 




•a 


l,n 




(1) 


l 




2,n 





(2) 



(2) 

2,71 



(4.17) 



so that both of them become nontrivial in the limit /x(= z 2 ) — > 1. Note that 
(1431) and (gSJ imply that 



(1) - z- x P n+1 *£ +1 + z- 1 (I - P n+1 Q n+1 ) $W 



2,n 
(2) 



■ 2,71+D 

■ 2,7i — * ~~ V - 1 n+l^l,n+l "T" * V- 1 ~~ - 1 n+lQn+i; ^2,n+l- 

Thus, the ratio between these quantities in the limit //(= z 2 ) — >■ 1 satisfies 



1 - 1) p n+ i* S+i + ^ (/ - p n+1 Q n+1 ) $f ) 



^2,71-1^2,77- 



/J— 5-1 
,(2)" 



(/ - p n g n ) |p n *g + (J - p n Q n ) $g 

-1 



(1) 



/i— >1 



(2) 

2,71 



P 1 P + $J, Z 

1 



,(1) 



GO -WD 

?i 2,n 



1,71 l,n 



/i— >i 



fj, — y 1 



Therefore, the formula for determining r n in (14. 16p can be replaced with 



-A+ (V 2) -W 1} 

^77. 1^2,71-1^2,71-1 



/J— >1 



(4.18) 



which uses the two linear eigenfunctions for (AL2). 
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We set 



(i) 



(i) 

2.n 



(ALl) 



(2) 



(2) 

2,71 



(ALl) 



(4.19) 



so that (cf. (Q7)» and (HTEjl ) 



In 
2,n 



(AL2) 



(2) 



(2) 
2.)i 



(AL2) 



(4.20) 



In view of ( 13.5)) . (14. 8 p and ( 14. 9p . this choice is indeed consistent with the 
boundary conditions (I4.15p . To be precise, we should use the explicitly 
time-dependent Jost solutions introduced in subsection 13. 4} which satisfy 
not only the Ablowitz-Ladik eigenvalue problem (13. ip but also the time- 
evolution equation ( 12. 2P ; however, this makes no difference in the limit // — > 1 
(cf. (13.31bp ). so the above choice is valid in formulas (I4.16P and ( I4.18p . Thus, 
with the aid of ( 13. lip , the solution of the space-discrete Kaup-Newell system 
(I4.14p can be expressed as 



k=0 



(ALl), 



n 



n + k) 



fc=0 



(ALl), 



n 



n + k) 



(4.21a) 



-a: 



n -1 



+E« 



(AL2) 
2 1 



n 



- l,n- 1 + k) 



k=0 



fc=0 



n 



l,n-l+k) 
(4.21b) 



Recall that the infinite sums in (14.21 p are assumed to be convergent; this 
is satisfied if the potentials in (ALl) and (AL2) decay sufficiently rapidly 
as n — > ±oo (cf. ( 13. lip ). To realize an exact linearization of the space- 
discrete Kaup-Newell system (I4.14p . we introduce new quantities /C(n, m) 
and /C(n, m) for m > n as 



/C(n, m) :- 
/C(n, m) :- 



I + J^K^ AL1 \n,n + k) 



k=0 



5>f AL1 W)> 



I + J2 Ki 2 AL2 \n-l,n-l + k) 

k=0 



Ki AL2 \n-l,s), 

s=m—l 
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so that q n = A+/C(n, n) and r n = A+/C(n, n). Let us derive the linear sum- 
mation equation for /C(n, m) from the Gel'fand-Levitan-Marchenko equa- 
tions f l3^2|) and f l3^9|) for (AL1). From f l^22|) . we have 

oo 

K^ AL1 \n,p) + J2 K { AL1 \ n ' n + j) FiAL1 \p + J + !) = O, p>n. 

j=0 

Thus, taking the sum with respect to p, we obtain the relation 



p=n+k 



(AL1), 



n 



p) = -E 



j=0 



E*! 



(ALl), 



n,s)- ir{ AL1) (ra ; -5) 

s=n+j+l 



x F( AL1 \p + j + l) 

p=n+k 



(4.22) 



From f)3.29p . we have 



K{ AL1) (n, s) = F^ AL1 \s) + Ki AL1) (n, n + k)F^ (s) 



k=0 



-E 

k=0 



E 

,p=n+k 



(ALl), 



n,p) 



F (ALl) (s + fc) _ F (ALl) (s + A , + 1) 



where (a variant of) the summation by parts formula is used. Thus, taking 
the sum with respect to s and using the fact that F^ Aljl \n) decays rapidly 
as n — >• +00, we obtain 



5>; ali W) 



k=0 



-E 



E ^. (al1> kp) 

fc=0 [_p=n+k 

Substituting (I4.22p into the last term and multiplying both sides from the 



£f<alD (s) 

s=m 

F {AL1 \m + k), m>n. 



left by 



'+£r=o^i (AL1 W+*o 



equation for /C(n,m): 
JC(n, m) 



we arrive at the linear summation 



00 00 



s=m j=0 k=0 



x F {AL1) (p + j + l)F (AL1) (m + k), m> 

p=n+k 



n. 



s > n, 
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which can be rewritten using (14. lip as 



oo oo 



/C(n, m) = F ( s ) + $3 53 n + A ~ ^( n ' n + 3 + ^ F ( n + J + k ) F ( m + *0> m ^ 

s=m j=0 fc=0 

In a similar way, we can derive the linear summation equation for K.(n, m) 
from the Gel'fand-Levitan-Marchenko equations (I3.22p and (I3.29p for (AL2) 
as 

OO M OO 

£(n,m) = f(AL2) ( s ) + EE[^ n+ ^-^ n ^' + 1 ) 

s=m— 1 j=0 k=0 

oo 

x J2 F^ AL2 \ P + j + l)F^ AL2 \m + k-l) 

p=n+k— 1 



n. 



OG 



OO oo 



F(s) - ^2Y1 [ lc ( n i n + 3) ~ K.{n,n + j + 1) F(n + j + k)F(m + k — 1), 

s=m— 1 j=0 fc=0 

m > n, 



where (14. lip is used. Combining the above results, we obtain a set of formulas 
for the solutions of the space-discrete Kaup-Newell system (14.141) . which tend 
to zero as n — > +oo, in the form [21J : 



q„ = A+JC(n,n), 



(4.23a) 
(4.23b) 



/C(n, m) 



K,(n, m) 



T{m) + ^ n+ 3)~ K ( n > n + 3 + !)] 

j=0 fc=0 

x [J"(n + j + k + 1) - T{n + j + k + 2)] T{m + k) - 7{m + k + 1) 

(4.23c) 



oo oo 



J*(m) - [^(n, n + j) - /C(n, n + j + 1) 

j=0 k=0 



m > n. 



x Tin + j + /c) — Tin + j + k + 1) [.F(m + fc) — .F(m + A; + 1)] , m>n. 

(4.23d) 

Here, the functions ^(n) and ^(n) are defined as Tin) := Yl^L n -i F ( s ) anc ^ 
•F(n) := Xl^=n respectively; they satisfy the same linear evolutionary 
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system as F(n) and F{n) (cf. (I3.38P ) and decay rapidly as n — > +00. Note 
that the set of formulas (I4.23P can provide the solutions for any flow of the 
space-discrete Kaup-Newell hierarchy if ^{n) and ^(n) satisfy the corre- 
sponding linear system. 

In the same way as for the other lattice systems, we can derive the mul- 
tisoliton solutions of the space-discrete Kaup-Newell system (14.141) from the 
set of formulas (14.231) ; this corresponds to the special case of reflectionless 
potentials for both (AL1) and (AL2). For simplicity, we assume that A(fi)^ 1 
and only have simple poles and set 



N N 

3=1 j=i 



(4.24) 



where the time dependences of Cj and Cj are given as (I3.34p and (13. 35 p . i.e. 

Cj{t) = Cj{0)^- 1}b+{1 -^ 1)a]t , Cj{t) = C i (0)e- [( ^- 1)6+(1 ~^ 1)a]t . 
We also set 

TV N 



K(n,m;t) = ^Gj(n,t)tMj 
3=1 



K(n, m;t) = J2 *)a#\ ( 4 - 25 ) 
3=1 



and substitute the expressions and (14^5) into ( I4.23cp and (I4.23d|) . 

Because \/jLj\ < 1 (j = 1, 2, . . . , N) and > 1 (j = 1, 2, . . . , N), we can 
evaluate the infinite sum to obtain linear algebraic systems for determining 
Qj and Tij, respectively. They can be written as 



fj n ft 



U 



IN 



u 



Nl 



u 



NN 



c 



N 



(4.26a) 



and 



V N1 



Vin 



[ Cifjii C2H2 ' ' ' CnUn ] 
(4.26b) 
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Here, all the entries in ( I4.26P are I x I matrices; the block matrices U = (Ujk)^ 
and V = (Vjfc)i<j,fc<Ar are defined as 

N 

Ujk := SjkJ 



^ (1 - Tir 1 ) (1 - uA (1 - 7J- 1 ) 



i=l 



1 



A** 



1 



M» 



and 



V 



w + e (1 "/ j)(1 ;y 1)(1 .r t) ™,wc tW . 



i=l 



1 _ 

Mi 



1 _ ^1 

Mi 



respectively. Thus, with the aid of (I4.23al) . (14.23bl) and (I4.25p . we obtain the 
multisoliton solutions of the space-discrete Kaup-Newell system (14.141) in the 
difference form: 



Qn(t) 




u 



IN 



u 



Nl 



u 



NN 



-1 














J 



(4.27a) 



r„(t) 





■ Vin ' 


-1 


" I ' 




_V m ■ 


■ ^NN 




I 


j 



(4.27b) 



When b = a*, the space-discrete Kaup-Newell system ( 14. 14j) allows not 
only the complex conjugation reduction r n = i&q*_u 2 but also the Hermitian 

conjugation reduction r n = io"g j |_ 1 y 2 , where a is a real constant [20]. Each 
reduction relates two originally uncoupled systems, (I4.14p with n G Z and 
(I4.14p with n G Z + 1/2, to provide a single equation with n G Z/2. Clearly, 
the value of a is nonessential, so we set o = 1 and consider the Hermitian 
conjugation reduction r n = i%_ x / T This reduction can be realized at the 
level of formulas (I4.23P by setting 

Hn) = i{F(n + \)}\ 

which is consistent with the time evolution (cf. (I3.38P with b = a*). In par- 
ticular, the iV-soliton solution of the space-discrete Kaup-Newell equation, 



Qn,t 



ail - \q n q\ A q n + a* (i + iq n -iql_i 



q n -i 



O, (4.28) 
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is obtained by setting b = a*, N = N and 

N = 4 Ci(t) = i fc W/iJ }\ j = 1, 2, .... AT 

in formula (I4.27aj) . where Cj(t) = C J .(0)e Kw_1)a * +(1 ~'*i" 1 > o ] t . The imaginary 
unit (roman i) should not be confused with the index of summation (italic 
i) in the definitions of Ujk and Vjk- The iV-soliton solution thus obtained 
is a space-discrete analog of the iV-soliton solution of the continuous matrix 
Kaup-Newell equation reported in [14] . Note that the square matrix equation 
(I4.28P can be further reduced to a vector equation by setting all but one of 
the columns (or rows) of q n to zero; its iV-soliton solution can be obtained 
in the same way as described in [5B] . 

5 Concluding remarks 

In this paper, we have developed the inverse scattering method associated 
with the matrix Ablowitz-Ladik eigenvalue problem and its applications to 
space-discrete analogs of derivative NLS systems. In particular, the most 
streamlined version of the inverse scattering method on a lattice is formu- 
lated, which can avoid redundant processes present in the existing literature. 
Thus, we are now able to understand the inverse scattering method for the 
Ablowitz-Ladik lattice as a direct discrete analog of the inverse scattering 
method for the continuous NLS system; in essence, the discrete case is no 
longer more complicated than the continuous case. Moreover, we can char- 
acterize the space-discrete derivative NLS systems using the potentials and 
linear eigenfunctions appearing in the Lax representation for the Ablowitz- 
Ladik lattice. On the basis of this characterization, we can solve the space- 
discrete derivative NLS systems by preparing two relevant copies of the in- 
verse scattering formulas for the Ablowitz-Ladik lattice and considering a 
Backlund-Darboux transformation between them. This provides a unifica- 
tion of the inverse scattering method for the NLS system and that for the 
derivative NLS systems in the discrete setting; such a unification is also pos- 
sible in the continuous case (see (EMU]). The multisoliton solutions of the 
space-discrete derivative NLS systems can be obtained in a straightforward 
manner within this unified framework; they reduce to the multisoliton so- 
lutions of the derivative NLS systems in the continuous space limit. Note 
that the space-discrete Kaup-Newell system allows the introduction of the 
potential variables with respect to the discrete spatial coordinate and can 
be rewritten locally in terms of these variables; our solution formulas reflect 
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this property accurately, that is, by construction any solution is written in 
the difference form using the forward difference operator. 

We assumed vanishing boundary conditions at spatial infinity, namely, as 
n — > +00 and n — > —00. In the case of matrix-valued dependent variables, 
this assumption imposes highly nontrivial conditions on the scattering data, 
which become almost trivial in the scalar case (cf. (I3.48P ). For simplicity, 
we derived such conditions in the reflectionless case of the potentials; how- 
ever, they are expected to be valid in the general case, because in the limits 
t — > ±00 the contribution of the continuous spectrum would become negligi- 
ble (cf. (I3.36P and (I3.37P with b = a*). Note also that our approach of solving 
the derivative NLS systems using the NLS eigenvalue problem is applicable 
under other boundary conditions that are amenable to the inverse scattering 
method or its generalizations. In addition, although we mainly considered 
the first nontrivial flows of the integrable hierarchies, our approach can be 
applied, with minor amendments, to the higher flows as well as the negative 
flows of the hierarchies. 

A Reduction to the vector modified Volterra 
lattice 

In this appendix, we consider the reduction of the matrix Ablowitz-Ladik 
lattice (12.31) to the vector modified Volterra lattice: 



Here, q n is a row vector of arbitrary dimension and ( ■ , ■ ) stands for the 
scalar product. The scalar (i.e., one- component) modified Volterra lattice 
was introduced by Hirota |56j and the two-component case was studied by 
Ablowitz and Ladik [6]. The vector generalization (lA.ip was recognized as an 
integrable system in the late 1990s (see, e.g., references in [ST]); it can also 
be considered as the simplest space-discrete analog of the vector modified 
KdV equation |32|l58]: 




As was shown in our previous paper [23] (also see [21]), the inverse scat- 
tering method formulated for the matrix Ablowitz-Ladik lattice (12. 3p can be 
applied to the vector modified Volterra lattice ( lA.ip by imposing suitable 
reduction conditions on the scattering data. Here, we revisit this problem 
and derive the reduction conditions in a more convincing manner. 




(A.l) 
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We introduce a set of 2 M 1 x 2 A/ 1 matrices {ei, e 2 , . . . , e 2 M-i} that are 
linearly independent and satisfy the anticommutation relations: 

{ej, e k } + := eje k + e k ej = -25 jk I. (A.3) 

Note that they are the generators of the Clifford algebra. We set the matrix- 
valued dependent variables Q n and R n as 

2A/-1 

Qn = q^I+ £ ti +l) e v (A.4a) 

3=1 
2M-1 

R n = - q Wl + tf + V (A.4b) 
Then, because of (IA.3I) . they satisfy the important relation: 

QnRn = RnQn = —(Qui Qn)I- 

Here, q n = (q n , . . . ,q n 2M ^)- Using this relation, it is easy to see that the 
choice (IA.4I) reduces the matrix Ablowitz-Ladik lattice (I2.3P with a = b = 1 
to the vector modified Volterra lattice flA.lj) . Lax representations involv- 



ing the generators of the Clifford algebra were introduced in the pioneering 
papers [SHEO], but our choice flA.41) is more efficient and useful because it 
contains the unit matrix /. As a natural analog of the complex conjugate, 
we define the Clifford conjugate for the linear span of {/, ei, e 2 , . . . , e 2 M-i} 
as 

2M-1 2M-1 

Qn = q^I - E Rn = " E ^ + %- 

3=1 3=1 

Note that Q n = Q n , R n = -Q n , Q n Q n = Q n Qn = (<?„, q n )I, etc. In short, 
the Clifford conjugate denoted by ^ changes the sign of the coefficients of 
{ei, e 2 , . . . , e 2 M_i}. This definition of the Clifford conjugate is very useful in 
the following discussion. 

Let us discuss how the reduction (1A.4j) constrains the scattering data for 



the matrix Ablowitz-Ladik lattice. The main role is played by the quantity 
P n = ^2,71^^ defined from the Lax representation, (12. ip and (12.21) . which 
satisfies the pair of discrete and continuous matrix Riccati equations (12.41) . 
We can show in an inductive manner that under the reduction (I A .40 . P n also 
takes its values in the linear span of {I, ei, e 2 , . . . , e 2 M-i}- We assume that 
the expression, 

2M-1 

Pn = - Pn 1) I+ Erf^^' 
3=1 
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is valid at some value of n and write its coefficients as p n := (ph , ■ ■ ■ ,Pn )• 
Then, noting the relations derivable from the anticommutation relations 

PnPn PnPn \PniPn)^- 1 

QnPn + PnQn = ~2{ Pn , q n )I, (A.6) 
QnPriQri (^QnPn PnQn) Qn PnQnQn 2 (p n , Qn)Qn (*?n' Qn)Pm 

(i ~ fiQuPn^j (I ~ pP n Qn) = (l + 2fi(p w <?„) + ^(Pn, Pn) (Qw Qn)) J ' 

we can rewrite (I2.4al) as 

P n+1 = (/ - nPnQnT 1 (jiPn + Qn] 



1 



(y^P-n Qn fJ> QnPnPn l^QnPnQi 



1 + 2p(p n ,q n ) + fi 2 (p n ,p n )(q 

ni Qn) 

= -i I o / rr 21 \l \ ^ ( X + («»> 9»» P « + I 1 + 2 ^(Pn> <7n> ~ ^ (Pm Pn)) Qr 

Thus, the expression (IA.5I) is also valid for P n+ i, and the coefficients satisfy 
the recursion relation written in the vector form: 

= M 1 + (QmQn))Pn ~ (1 + 2 MPn ? gn) ~ ^ (Pn. Pn) ) <?n 
^ 1 + 2/i(p n ,q n ) + H 2 (Pn,Pn)(Qn,Qn) 

In a similar way, we can show using fl2.4a|) that if the expression (1A.5|) is valid 
for P n+ \, then it is also valid for P n . Therefore, under a suitable boundary 
condition on P n , such as lim^-oo P n = O or \im n ^, +00 P n = O, so that the 
expression AA.5D is valid at the boundary, P n indeed takes its values in the 
linear span of {J, ei, e2, . . . , e2M-i} f° r & U n E Z. Moreover, using (12.4bl) . we 
can easily check that this property is preserved under the time evolution. 
Substituting the expressions 

Pn Q"n,Pn b n Pn+li 

Qn Rn &nPn &nPn+l 

into (12.4aj) . we can express the scalar coefficients a n and b n in terms of p n and 
p n+ i- Thus, q n can be written explicitly using p n and p n+1 , which defines a 
Miura map to the vector modified Volterra lattice flA.10 . The corresponding 
modified system, i.e., the closed differential- difference equation for p n can be 
obtained from (I2.4bl) . but we do not present it here (cf. [571I6T] ). 



49 



Because (|2.4a[) can be rewritten as 



it is possible that fi~ n P n for — 1 has finite limits as n — > ±00 under the 
decaying boundary conditions on Q n and R n (cf. (13 .2p ). In particular, the 
above difference equation with the boundary condition lim^^oc fi~ n P n = O 



defines the other boundary value lim r 



/i n P n , which belongs to the linear 



span of {/, ei, e2, . . . , e2j\/-i}. The defining relation (I3.7ap of the scattering 
data together with (13. 5p implies that on the unit circle = 1, 



z n I O 
O z~ n I 



I O 

O fx~ n I 

I 

o 



z n ip n A + 



fi n I O 
O I 



as n —> —00, 



as n — > +00. 



z~ n ip n B 



Thus, if we set fi~ n P n as fi~ n P n = (z~ n (p2,n)(z n (pi tn )~ 1 using the above linear 
eigenfunction, we find that lim n _j. +00 {i~ n P n = takes its values 

in the linear span of {/, e±, . . . , e2M-i}- More precisely, we should use 
the explicitly time-dependent Jost solution <pn as introduced in f)3.31ap . but 
this makes no difference in the above discussion. In addition, the established 
property of -B(/i)A(/i) _1 is preserved under the time evolution. 

For < 1, we can still consider the difference equation f)2.4al) for P n with 
the boundary condition linin^-oo P n = O by setting P n = (z n (p2,n)(z n (fii^)- 1 ; 
however, in this case the other boundary value lim n _> +00 \i~ n P n does not exist 
in general. Thus, we need to take a more delicate limit for n and /x to extract 
meaningful information from P n . Let us evaluate the n — > +00 behavior of 
the Jost solution [z n (j) n \{ii) in the neighborhood of n = jij < 1) at which 
t4(/z) _1 has a simple pole. Recalling the results in subsections 13.21 and 13.31 
we obtain that 



[z n n ](/i) x - iij)A(fi) 



-1 



- H) {I + o(l)} + o(j$) 
(ji-tij) o(l) + C jf x] + o(^) 



Thus, if we take the limit of n — > +00 and — > fij while maintaining the 
balance condition |/i — ~ |//"| , we arrive at the desired formula: 



/i — fij 



C 3l x] + o(l). 
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Hence, Cj also takes its values in the linear span of {/, ei, e 2 , . . . , e 2 A/-i}; 
naturally, this property is preserved under the time evolution. 

Combining the above results, we conclude that F(n,t) as given in (I3.36P 
with a = b = 1 (or its generalization corresponding to the case of higher order 
poles of A(fi)^ 1 ) can be expressed in the form: 

2M-1 

F(n,t) = f®(n,t)I- (A.7) 

0=1 

We remark that the quantity V n := $i in \l/^ defined from another linear 
eigenfunction of the Ablowitz-Ladik eigenvalue problem (12. II) satisfies 

Q n = l±~ l V n - V n+ i + [i~ 1 V n R n V n+1 , 

or equivalently, 

V n+ i = (I- ^VnRnY 1 {yT X V n - Q n ) . 

Because Q n and R n in (1A.4|) are related as R n = —Q n , we can identify V n 
with the Clifford conjugate of P n as 

Here, we assume that the boundary conditions for V n and P n are compatible 
in the above identification, e.g., \im n ^ ) ._ 00 V n = lim^.oo P n = O. Thus, if 
wc set "T^fi — (z n <f) ljn )(z n <p2,n) 1 using the Jost solution <p n as introduced in 
(I3.5ap (or the explicitly time-dependent one <f>n in (I3.31a|) ) and compare it 
with P n = (z n (j)2 } n){z n 4'i,n)~ 1 , we obtain the following relations: 



B(n)A(fi) 1 on \fi\ = 1 is the Clifford conjugate of B(fj, 1 )A( / u 



i\— i 



• ix j — — , Cj = ^Cj, j = 1,2, ... , N(— N), up to renumbering. 

Therefore, F(n, t) as given in (I3.37P with a = b = 1 is equal to the Clifford 
conjugate of F(n,t): 

F(n,t) = F(n,t) 

2M-1 

= fV(n,t)I+ £ f U+1 \n,t)e 3 . (A.8) 
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Here, the vector f n (t) := (f^(n, t), . . . , f( 2M \n, £)) satisfies the linear evolu- 
tion equation f n t = f n+1 — (cf. (13.381) ) and decays rapidly as n — > +00. 

In fact, (1A.7I) and (1A.8j) provide not only necessary but also sufficient 
conditions for the corresponding potentials Q n and R n to be expressed as 
(1A.4I) without using e^e^, eitjtk-, etc. 



Proposition A.l. If F(n,t) and F(n,t) are given as (IA.7j) and flA.81) . then 



the potentials Q n and i? n reconstructed from the set of exact linearization 
formulas ( I3.39P can be expressed in the form (IA.4p . 



Before proving Proposition \A.1\ we need to state one proposition and 
two lemmas. We first show that the set of formulas (I3.39P realizes an exact 
linearization of the matrix Ablowitz-Ladik lattice. 

Proposition A. 2. Suppose that F(n, t) and F(n, t) satisfy the pair of linear 
evolution equations ( I3.38P and decay sufficiently rapidly as n — > +00. Then, 
the Liouville-Neumann-type series for the potentials Q n and R n defined by 
formulas (I3.39P solve the matrix Ablowitz-Ladik lattice (I2.3P exactly. 

Remark. Only the case a = b = 1 is relevant to the vector modified Volterra 
lattice (lA.ip , but we find it more convenient to prove Proposition IA.2I for 
general values of a and b. 

Proof of Proposition [X7B . Using (I3.39p . we obtain the Liouville-Neumann- 
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type series for Q n and R n in the form: 

oo 

Q n = F{n)- F{n + h)F{n + i 1 +i 2 + l)F(n + i 2 + l) 

ii,i2=0 

oo 

+ F(n + i 1 )F(n + ix + i 2 + l)F(n + i 2 + i 3 + l)F(n + i 3 + i 4 + l)F(n + u + l) 



fc=0 

oo 

R n = -F{n) + F(n + i x )F{n + h + i 2 + l)F(n + i 2 + 1) 



h,i 2 =0 



- Y F{n + i x )F{n + ii+i 2 + l)F(n + i 2 + i 3 + l)F{n + i 3 + u + l)F(n + i A + l) 

*l)*2,«3i*4=0 
_|_ ... 

OO 

= ^(-l) fc+1 7^. (A.9b) 

k=0 

Here, 

:= F(n), Tlf := F(n), (A.lOa) 
and Ql fc) and TZ^ for fc > 1 are defined as 

oo 

2i fc) : = Yl F(n + i 1 )F(n + i 1 +i 2 + l)F(n + i 2 + i 3 + l)F(n + i 3 +i4 + l) 

ii,i2,— ,«2fc=0 

x • • • F(n + i 2 fe-2 + «2fc-i + + «2fc-i + i2k + l)F(n + i 2fc + 1), 

(A.lOb) 

oo 

K n k) := F(n + i 1 )F{n + i l + i 2 + l)F(n + i 2 + i 3 + l)F{n + i 3 + H + l) 

ii,i2,...,*2fc=0 

x • • • F(n + i 2k - 2 + z 2 fc-i + l)F{n + i 2k _ x + 2 2fe + l)F(n + i 2k + 1). 

(A.lOc) 

Note that Q n k) and 1Z n k) are "polynomials" of degree 2k + 1 in F and F with 
their arguments bounded below but unbounded above. We only consider the 
region of (n, t) where the series in (IA.9j) are absolutely convergent and admit 
termwise differentiation by t. 
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In the following, we use the shift operator E n as well as its inverse E n 1 
defined as 

E n f{n) := f(n + 1), E?f(n) := f(n - 1), 
and the forward difference operator in each index of summation i a , 

a;„. /•(..../ ) :=/(..., i a + i,. ..)-/(•• .,*«,...)■ 

Then, using (I3.38p . we can compute the time derivative of for > 1 to 
obtain 

[<9 t - aE n + + (a - b)l] Q« 
= a J F(n + ij + l)F(n + ^ + i 2 + l)F(n + 1 2 + z 3 + 1) • • • F(n + i 2fc + 1) 

ii,i2,.--,«2fc=0 

- F(n + h)F{n H + ^)^(« + i 2 + i 3 + 1) • ■ • F(n + i 2k + 1) 
+ ... 

- F(n + i^Fin + i 1 +i 2 + !)■■■ F (ji±j^k-A±i2k)F(n + i 2k + 1) 

+ F(n + h)F(n + i x + i 2 + 1) • ■ • F(n + i 2fe _i + i 2fe + l)F(n + i 2 fc + 2) 

- F(n + h± D F (n + *i + *2 + 2) ■ • • F(n + i 2k -i + kk + 2)F(n + i 2k + 2) J 

CO f 

+ b <-F(n + i 1 -l)F{n + i 1 +i 2 + l)F{n + i 2 + i 3 + l)---F{n + i 2k + l) 

ii,iai—,»afc=0 

+ F(n + h)F(n + h + i 2 + 2)F(n + i 2 + £ 3 + 1) ■ • • F(n + i 2fc + 1) 

- F(n + i x )F{n + i x + i 2 + l)F(n + i 2 + ja) • • • F{n + kk + 1) 
+ ... 

+ F(n + h)F(n + i x + i 2 + 1) • • • ^(j^tiB^l±i3i±S^( n + ^ 2fc + X ) 

- F(n + i x )F{n + n + % 2 + 1) • • • F(n + i 2k ^ + i 2k + l)F(n+j2k) 

where the arguments shifted from the original ones in Q n are underscored 
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with a wavy line. We can further rewrite it as 



a 



*l:*2v.*2fc=0 



+ A 



13 



F(n + h)F(n + + i 2 )F(n + i 2 + i 3 + 1) • • • F(n + i 2k + 1) 



F(n + ii)F(n + ii + i 2 + l)F(n + i 2 + i 3 + l)F(n + i 3 + u) • • • F(n + i 2fe + 1) 



+ 



»2fc-l 



F(n + ii) ■ ■ ■ F(n + i 2k _ 2 + kk-i + l ) F (jt+j3.t±+J2t) F ( n + i2k + ^ 
A+ a _ x |>(n + h)F{n + h + i 2 + 1) • • • F(n + i 2a _ 2 + i 2a ^ + 1) 

x F{n + i 2a -i + i 2a + l)F(n + i 2a + Wi + 2) • • • F(n + i 2k + 2) 



a=l 



U,»2,— ,»2k=0 



A+ 



F(n + ii)F(n + ii + i 2 + l)F(n + i 2 + i 3 ) • • • F(n + i 2k + 1) 



+ ■■■ 

+ a: 



F(ra + ii) ■ ■ ■ F(n + h + U + l).F(ra + i 4 + i 5 ) • • • F(n + i 2k + 1) 



*2fc 



k 

0=1 



F(n + ii)F(n + ii + i 2 + 1) • • • F(n + i 2k -i + %i k + VH^+jg) 

Fin ■ /'; 1 )/'( /' • %\ + 12) ■ ■ ■ F{n + « 2 /3-i + i2p)F{n + i 2( g + «2/3+i) 
x F(n + i 2/ 3+i + Vf 2 + !)••• F(n + i 2k + 1) } 



= a S A £«-i A £> [-^ + «l)^ + «l+«2 + l)---i 7, (n + i 2a -2+i2a-l + l) 

11,12,— >»2fc=0 !<"</3<fc 

X ^(^J^a-H^j^ + ^2a + Wl + !)••• ^(n + 12/3-1 + «2/3 + 1) 

x F(n + % 2f} + i 2 /3+i + ^^ftij^+H^ • • • ^(^ti25jil2} 

oo 

ii,i2,---,*2fc=0 l<o</3<fe L 

X ^(^J^J^a;^^^" + «2a + Wl + !)••• ^(n + 120-1 + ^2/3 + 1) 

x F(n + i 2/3 + i 2/3+ i)F(n + i 2 p +1 + 22/3+2 + !)••• F(n + i 2k + 1) 



which is equal to 



a E - E F(n + i 1 )F(n + i 1 +i 2 + l)---F(n + i 2a _ 2 + l) 

l<a<f)<k U,...,i2a-2=0 

oo 

X E F ( n + i2a)F{n + l2a + *2 a +l + 1) • • • F(n + i 2 /3-l + 1) 

'2a, ...,12/3-1=0 

OO 

x F(n + t w+1 + l)F(n + i 2f3+1 + i w+2 + 2)---F(n + i 2k + 2) 

»30+i>— i«2fe=0 



+ & E E ^(n + H - + ^ + z 2 ) • • • F(n + i 2a _ 2 ) 

l<a</3<fc ii,...,i2 a -2=0 
oo 

x E F(n + i 2 a)F(n + i 2a + i 2a+1 + 1) • • • F(n + i 2 p-i + 1) 

«2a, ■••1*2/3-1=0 
oo 

x E F(n + z 2j9+1 )F(n + z 2/3+ i + ^+2 + 1) • • • F(n + i 2k + 1) 

*2/3 + l,-..,*2fc=0 



^<a<P<k l<a</3<k 



Thus, we obtain 

[d t - aE n + bE- 1 + (a- b)l] 

= -a E eSr^^oaf +6 E sfciM^Q^, *>i. 

(A.lla) 

Similarly, we also obtain 

[dt - bE n + aE- 1 + (b- a) I] 

= -b E ^Qt^^ + a E ^Qr^^), fc>l. 

l<a</3<fc l<a</3<fc 

(A.llb) 



Because of flA.9j) . flA.lOaj) and (ET3H1) . relations (lA.lljl imply that 

Qn,t ~ Q-Qn+i + bQ n -i + (a — 6)Q n = —aQ n R n Q n+ i + bQ n -\R n Q T , 
Ft n ,t — bR n+ i + ai?„_i + (6 — a)i? n = —bR n Q n R n+ i + aR n _iQ n R n 
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These are exactly the equations of motion for the matrix Ablowitz-Ladik 
lattice (J23D- □ 

blowing, we set a = b = 1, which simplifies (13.381) to 

(A.12a) 
(A. 12b) 



In the following, we set a = b = 1, which simplifies (13.381) to 
d t F(n, t) = F(n + 1, t) - F(n - 1, t), 
dtF(n, t) = F(n + 1, t) - F(n - 1, t). 



Here, F(n,t) and F(n,t) are required to decay rapidly as n — > +oo. For 
the moment, we do not impose the conditions (1A.7j) and (1A.8j) . and consider 
the general case where F(n,t) and F(n,t) are independent matrix functions. 
Then, we can rewrite flA.llj) as recurrence relations for Qn and TZ^ that 
were originally defined as flA.lOp . i.e. 



QW = (J + E n d t - E 2 n ) 



2N-1 



K^ = (l + E n d t -Eiy l 



(0) 



^n+l 'Si+1 ~r / „ 'Si+1 

l<a</3<fc !<«</3<fc 

(A.13a) 



(q-1) n (/3-a)-r> (fc-/3) 
'Si 

l<a</3<fe 

(A.13b) 



„ - F(n,t) and 7?-™ = F(n,t). Here, the inverse operator 



for > 1, and Q 
(J + - ££)~ A is defined using the Maclaurin series as 

oo 

(/ + E n d t - E 2 n ) 1 :=J2K (E n - d t y , 

j=0 



where the action of dt on F(n,t) and F(n,t) is given by (1A.12j) . 

In fact, the linear operator / + E„d t — E% has a nontrivial kernel that 
contains F(n, t) and F(n, t) as well as their spatial/temporal shifts, so it is not 
invertible in general; however, this does not cause any problem in obtaining 
(1A. 13|) . To confirm this, it is sufficient to note that no "polynomials" of 
degree 2k + 1 (k > 1) in F and F can vanish by the action of I + E n d t — E^, 
as long as their spatial arguments are bounded below. 

Lemma A. 3. For general F(n) and F(n) satisfying the linear evolution 
equations (IA. 12j) . assume that the following equality is valid: 

oo 

(I + E n d t -El) E li 1 i 2 ...i 2k+1 F(n + i 1 + a 1 )F(n + i2 + a 2 )F(n + i 3 + a 3 ) 

*l>«2,--->«2fc + l=0 

x • ■ ■ F(n + i 2k + a 2k )F(n + i 2k +i + a 2k+ i) = O. 
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Here, 7i 1 j 2 ...i 2fc+1 E C and a,- E Z> are constants, and k > 1. Then, the equal- 
ity must be trivial, i.e. 

7nta-»2fc+i = 0' ^1)^2, • • • , ^2fc+i € Z> . 



This lemma can be easily verified by ordering the nonzero terms involved 
in the summation, e.g., using sums of the arguments such as 

(n + i x + «i) + (n + z 2 + a 2 ) + (n + z 3 + a 3 ) H h (n + i 2 fc+i + «2fe+i), 

(n + z 2 + a 2 ) + (n + i 3 + a 3 ) H h (n + z 2 fc+i + "2/t+i), etc. 

Note that both -E n <9 4 and i?^ increase the values of the arguments in each 
term, while / leaves them invariant. Thus, the equality implies that the coef- 
ficient of the first term in the ordering, which has the minimum values of the 
arguments, must be zero. Therefore, there is no first term and all 7i 1 i 2 .„i 2fe+1 
must vanish. □ 

Lemma DOl guarantees that the difference between Qn defined as flA.lObl) 
and the right-hand side of flA.13al) does not belong to the kernel of I + E n dt — E\. 
Thus, (IA.13al) is indeed an exact equality; the same applies to (IA.13bl) . 

On the basis of the recurrence relations (IA.13I) . we can prove Proposi- 
tion IA.1I by induction. To this end, we need to use one lemma, which is a 
direct consequence of the anticommutation relations (1A.3j) for the generators 
of the Clifford algebra (cf. (tPjl ). 



Lemma A. 4. If the square matrices X and Y take their values in the linear 
span of {/, d, e 2 , . . . , e 2A /-i}, then 

XY + YX 

is a scalar matrix and coincides with XY + YX. Thus, up to an overall 
factor, it can be considered as the definition of an inner product. Here, 
denotes the Clifford conjugate, which changes the sign of the coefficients of 
{ei,e 2 , . . . ,e 2A /-i}. 



Proof of Proposition Because of (1A.9[) . we need only show that, for all 
k E Z> , 

Qn^ and IZn take their values in the linear span of 

{/, ei, e 2 , . . . , e 2 A/_i} and satisfy the Clifford conju- (♦) 

gation relation = Qn\ 
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Note that this is true for k = (cf. (TjOJ, ([Q|) and flA.lOaj) ). We assume 
that (4k) is true for k < L and then proceed to show (4k) for k = L + 1. 
Lemma IA.4I implies that the quantity appearing in ( 1A.13aj) with k = L + 1, 



/3 
a=l 



n+1 /v n+l 



^(a-lJ-nOS-a) , r) (^-a) 7? (a-i; 

2 2-~i L n+1 "+ 1 ^n+1 'Si+1 
a=l 

1 ^ 

2 [^n+l ^n+1 ^n+1 %i+l 



1 </3 < L + 1, 



a=l 



is a scalar matrix and coincides with the quantity appearing in ( lA.13bj) with 
k — L + 1, 



E (a-l) n (P-a) _ 1 I" ( a -l) n (/3-a) (/3-a) n (a-l) 

'Si+1 *in+l — o / v 'Si+l *Wl ~r 'Si+1 ^n+1 



a=l 



a=l 

/8 



El A(a-1) n C8-a) , fjfp-a) n (a- 
^n+1 ^n+1 i" ^n+1 ^n+l 



08-0)^,(0-1) 



1</3<L + 1. 



a=l 



Similarly, we also have 



L+1 



Er^ip-a) n {L+X-p) 
'Si+1 ^n+1 



L+1 

£s<' 

/3=a 



+1 'Si+1 



scalar, 



1 < a < L + 1. 



Thus, the recurrence relations (1A.13|) imply that (4k) holds true for k = L + 1. 
Therefore, by complete induction, (4k) is true for all k G Z> . □ 



Remark. By considering the continuous space limit of Proposition IA.14 we 
can obtain similar results on the corresponding reductions of the continuous 
matrix NLS hierarchy. In particular, we can identify the reduction condi- 
tions on the scattering data to solve the vector modified KdV equation (|A.2[) 
and the vector sine-Gordon equation |59j by the inverse scattering method 
(see |62j). In addition, it is an easy task to consider a continuous analog of 
Proposition lA.2l that realizes an exact linearization of the matrix NLS system. 

Proposition IA. II demonstrates that under the restrictions OA .70 and (1A.8|) 
on F(n,t) and F(n,t), we can solve the vector modified Volterra lattice 
(lA.ip by the inverse scattering method. Moreover, if the vector dependent 
variable q n is real-valued, we need to restrict the components of the vector 
f n (t) = (/W(n,t), . . . , f( 2M 1(n,t)) appearing in fjA~7]) and CO]) to be real- 
valued. Recall that in the reflectionless case of the potentials, if A(/i) _1 only 
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has isolated simple poles, F(n,t) is given as (cf. (I3.36P with a = b = 1) 



N 

F(n,t) = -J2CM^ {H ~" ;1)t , (A.14) 

3=1 

where < \fij\ < 1. Thus, the configuration of fi 2 , . . . , (jln} must be sym- 
metric with respect to the real /i-axis, that is, they either take real values or 
occur in complex conjugate pairs. Up to a reordering, they can be classified 
into the following three types: 

(I) ^2k-i = P* 2k = a k e w \ < a k < 1, < 9 k < tt, k = 1, 2, . . . , N 1: 

(II) < fij < 1, j = 2Ni + 1, . . . , 2Nt + N 2 , 

(III) -1 < fij < 0, j = 2N 1 +N 2 + l,..., 2N t + N 2 + N 3 (= N). 

Each of the corresponding matrices {Ci(0), 6*2(0), Cat (0)} takes its values 
in the linear span of {/, ei, e2, . . . , e 2 M-i}, wherein the coefficients are all 
real for type (II) and type (III). For type (I), the coefficients associated with 
fi2k-i and those associated with fj, 2k form a complex conjugate pair so that 
the coefficients of {I, ei, e 2 , . . . , e 2 M-i} in HA. 141) become real- valued [231121] . 

The nature of the three types (I) — (III) can be unveiled by constructing 
the corresponding one-soliton solutions. Type (II) provides the trivial vector 
analog of the one-soliton solution of the scalar modified Volterra lattice, 

sinha 

0n W = — n n , ■ 1 — r, ft u ) {u,u) = l. 

nW cosh[an + 2(sinha)t + 5] 

Type (III) provides a very similar solution, 

/ \ ( — 1)™ s i nn a 1 \ 

Quit) = — T-f of ■ r rr— »■ («, «) = 1, 

cosh [an — 2(smha)t + d\ 



which reflects the form-invariance of the vector modified Volterra lattice f lA.ip 



under the transformation q n — >■ (— l) n q n , t — > —t. For type (I), we need to 
impose an additional condition to exclude a breather solution so that a pure 
soliton solution with a time-independent profile can be obtained 
Thus, the one-soliton solution of type (I) is given by 



Qn(t) 



£gi/3n+2i(cosh a sin 0)t _j_ q*q~ i/3n— 2i(cosh a sin j3)t 

cosh [an + 2(sinh a cos 0)t + 5] 



with (c, c) = and 2(c, c*) = (sinha) 2 . This is indeed the most general one- 
soliton solution, because it reduces to type (II) and type (III) in the limit 
f3 — > and — > tt, respectively. 



60 



Proposition IA. II enables us to formulate the inverse scattering method for 
the vector modified Volterra lattice flA.ip . as well as the relevant continuous 
systems, with satisfactory rigor. As by-products, it can also generate nu- 
merous nontrivial identities for some rational functions of {/ii,/i2, . . . ,^n}- 
Indeed, for F(n,t) and F(n,t) given as (IA.7j) and (1A.8j) . the coefficients of 



e j e k (j 7^ k), CiCjCk (i ^ j 7^ k 7^ i), etc. in Q« for k > 1 defined as flA.lObj) 



must vanish (cf. (Jfr)). For the reflect ionless case of the potentials, we can 
express these coefficients explicitly in terms of the scattering data, which 
indeed provide nontrivial identities. If only has isolated simple poles, 

F{n) and F(n) are given as 

N N 

F( n ) = -Y,Cjti, F(n) = -Y,CjHl (A.15) 

3=1 3=1 

Here, {/ii, fi 2 , ■ ■ ■ ,^n} are pairwise distinct and satisfy < \fij\ < 1, and ~ 
denotes the Clifford conjugate. The time dependence of Cj and Cj is irrel- 
evant in this context and thus is omitted. Substituting (1A.15I) into (lA.lObj) 



and computing the multiple sum with respect to z 1( i 2 , ■ ■ ■ , i 2 k, we obtain 

N n n i i n+l i i n+l i i n+l ■ ■ ■ ii n+1 n n+1 r r r r ...r r 

I'.i.l'.i: f'h /'./: /'/,/ /'/,/ • ( J A 32 l ;.: ( 34 L '32k L '32k+l 

3i>h,—>hh+i=l 



Qik) = _ j2 — 



for k > 1. Thus, we can state the following proposition. 

Proposition A. 5. For {ei, e 2 , ■ ■ ■ , e2M-i} satisfying the anticommutation 
relations eje^ + e^ej = —25jkl, all the coefficients of quadratic or higher 
terms such as (j 7^ k) and titjtk (i 7^ j 7^ k 7^ i) in the quantity, 

r** r** 

ST^ il ( - y J2 ( - y i3°j4 ' ' ' ( - y j2fc°i2fc+l 

r, , , ; A*il (1 - /^lMj 2 ) (1 - Mi 2 ^ 3 ) (1 - ^3^4) •••(!- ^32k^32k+l) ' 

Ul J2,-J2tfl) 
= {jl,l2,.-->*2fc+l} 

vanish identically. Here, {zi, z 2 , . . . , i 2 k+i} is any (repeated) combination of 
positive integers (1 < i\ < i 2 < ■ ■ • < i 2 k+i) and 



2M-1 2M-1 



(a+l) 



Because {cf \ c?\ . . . , cj 2M ^} for j e {z 1; z 2 , . . . , «2fc+i} at a fixed time can 
be chosen arbitrarily, each identity obtained this way splits into as many 
identities as the number of different products of {c^ } involved in the identity. 
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For example, when i 2 , . . . , 22/c+i} ar e pairwise distinct, we can simply 
set them as {1, 2, . . . , 2k + 1}. Thus, the coefficients of the highest products 
of {ej} such as eie2 • • • e2k+i provide the identity: 



sgn 



1 2 ... 2k + l 
ji h ■■■ hk+i 



|J1 J2v J2fc+l} 
= {l,2,...,2fc+l} 

Considering the lower products such as eie2 • • • eik-, • • • ^-i, • • ., eie2, 
we obtain various extensions of the above identity, wherein the sign of the 
permutation is replaced with more elaborate functions that still take values 
in {+1,-1}. 

Proposition IA.5I can be generalized in many different directions. First, 
we can consider the case where v4(yu) _1 also has second or higher order poles. 
Thus, pJf in flA.15j) can be replaced with more general functions decaying 



as n — > +oo. Second, instead of using Proposition IA.1I designed for solving 
the vector modified Volterra lattice flA.ip . we can consider corresponding 
results for other integrable systems having a Lax representation of the same 
type. That is, any integrable system with vector dependent variables can 
provide a similar result if it can be obtained as a reduction of a matrix- valued 
integrable system using the generators of the Clifford algebra. In particular, 
starting with the solution formulas for the matrix generalizations of the NLS, 
derivative NLS, modified KdV as well as their discrete analogs [21], we can 
obtain various interesting variants of Proposition IA.5I 



References 

[1] V. E. Zakharov and A. B. Shabat: Exact theory of two-dimensional 
self-focusing and one- dimensional self-modulation of waves in nonlinear 
media, Sov. Phys.-JETP 34 (1972) 62-69. 

[2] V. E. Zakharov and A. B. Shabat: Interaction between solitons in a 
stable medium, Sov. Phys.-JETP 37 (1973) 823-828. 

[3] M. J. Ablowitz, D. J. Kaup, A. C. Newell and H. Segur: Nonlinear- 
evolution equations of physical significance, Phys. Rev. Lett. 31 (1973) 
125-127. 

[4] M. J. Ablowitz, D. J. Kaup, A. C. Newell and H. Segur: The inverse scat- 
tering transform — Fourier analysis for nonlinear problems, Stud. Appl. 
Math. 53 (1974) 249-315. 



62 



T. Tsuchida: A systematic method for constructing time discretizations 
of integrable lattice systems: local equations of motion, J. Phys. A: Math. 
Theor. 43 (2010) 415202. 

M. J. Ablowitz and J. F. Ladik: Nonlinear differential-difference equa- 
tions and Fourier analysis, J. Math. Phys. 17 (1976) 1011-1018. 

V. S. Gerdzhikov and M. I. Ivanov: Hamiltonian structure of multi- 
component nonlinear Schrddinger equations in difference form, Theor. 
Math. Phys. 52 (1982) 676-685. 

D. J. Kaup and A. C. Newell: An exact solution for a derivative non- 
linear Schrddinger equation, J. Math. Phys. 19 (1978) 798-801. 

H. H. Chen, Y. C. Lee and C. S. Liu: Integrability of nonlinear Hamilto- 
nian systems by inverse scattering method, Phys. Scr. 20 (1979) 490-492. 

M. J. Ablowitz, A. Ramani and H. Segur: A connection between non- 
linear evolution equations and ordinary differential equations of P-type. 
II, J. Math. Phys. 21 (1980) 1006-1015. 

V. S. Gerdjikov and M. I. Ivanov: The quadratic bundle of general form 
and the nonlinear evolution equations. II. Hierarchies of Hamiltonian 
structures, Bulg. J. Phys. 10 (1983) 130-143 [in Russian]. 

A. Kundu: Landau-Lifshitz and higher-order nonlinear systems gauge 
generated from nonlinear Schrddinger-type equations, J. Math. Phys. 25 
(1984) 3433-3438. 

T. Tsuchida: Soliton solutions to integrable lattices of the derivative NLS 
type (2007), available at |http:/ /poisson.ms.u-tokyo.ac.jp/~tsuchida/| . 

T. Tsuchida: Exact solutions of the potential deriva- 
tive nonlinear Schrddinger equation (2008), available at 



http://poisson.ms.u-tokyo.ac.jp/ tsuchida/ . 



M. J. Ablowitz: Lectures on the inverse scattering transform, Stud. 
Appl. Math. 58 (1978) 17-94. 

M. J. Ablowitz and H. Segur: Solitons and the Inverse Scattering Trans- 
form (SIAM, Philadelphia, 1981). 

T. Tsuchida: Systematic method of generating new integrable systems 
via inverse Miura maps, J. Math. Phys. 52 (2011) 053503. 



63 



[18] R. M. Miura: Korteweg-de Vries equation and generalizations. I. A 
remarkable explicit nonlinear transformation, J. Math. Phys. 9 (1968) 
1202-1204. 

[19] P. D. Lax: Integrals of nonlinear equations of evolution and solitary 
waves, Commun. Pure Appl. Math. 21 (1968) 467-490. 

[20] T. Tsuchida: Integrable discretizations of derivative nonlinear 
Schrddinger equations, J. Phys. A: Math. Gen. 35 (2002) 7827-7847. 

[21] T. Tsuchida: New reductions of integrable matrix partial differential 
equations: Sp{m) -invariant systems, J. Math. Phys. 51 (2010) 053511. 

[22] V. E. Zakharov and A. B. Shabat: A scheme for integrating the nonlinear 
equations of mathematical physics by the method of the inverse scattering 
problem. I, Funct. Anal. Appl. 8 (1974) 226-235. 

[23] T. Tsuchida, H. Ujino and M. Wadati: Integrable semi-discretization of 
the coupled modified KdV equations, J. Math. Phys. 39 (1998) 4785- 
4813. 

[24] T. Tsuchida, H. Ujino and M. Wadati: Integrable semi- discretization of 
the coupled nonlinear Schrddinger equations, J. Phys. A: Math. Gen. 32 
(1999) 2239-2262. 

[25] M. J. Ablowitz and J. F. Ladik: Nonlinear differential-difference equa- 
tions, J. Math. Phys. 16 (1975) 598-603. 

[26] F. Calogero and A. Degasperis: Nonlinear evolution equations solvable 
by the inverse spectral transform. I, Nuovo Cimento B 32 (1976) 201- 
242. 

[27] O. O. Vakhnenko: Solitons in parametrically driven discrete nonlinear 
Schrddinger systems with the exploding range of intersite interactions, 
J. Math. Phys. 43 (2002) 2587-2605. 

[28] S.-C. Chiu and J. F. Ladik: Generating exactly soluble nonlinear dis- 
crete evolution equations by a generalized Wronskian technique, J. Math. 
Phys. 18 (1977) 690-700. 

[29] F. Kako and N. Mugibayashi: Complete integrability of general nonlinear 
differential- difference equations solvable by the inverse method. II, Prog. 
Theor. Phys. 61 (1979) 776-790. 



64 



[30] R. I. Yamilov: Classification of discrete evolution equations, Usp. Mat. 
Nauk 38:6 (1983) 155-156 [in Russian]. 

[31] V. E. Adler, S. I. Svinolupov and R. I. Yamilov: Multi- component 
Volterra and Toda type integrable equations, Phys. Lett. A 254 (1999) 
24-36. 

[32] S. I. Svinolupov and V. V. Sokolov: Vector-matrix generalizations of 
classical integrable equations, Theor. Math. Phys. 100 (1994) 959-962. 

[33] A. V. Mikhailov and A. B. Shabat: Integrability conditions for systems 
of two equations of the form u t = A{u)u xx + F{u, u x ). II, Theor. Math. 
Phys. 66 (1986) 31-44. 

[34] A. V. Mikhailov, A. B. Shabat and R. I. Yamilov: The symmetry ap- 
proach to the classification of non-linear equations. Complete lists of 
integrable systems, Russian Math. Surveys 42:4 (1987) 1-63. 

[35] Y. Ishimori: An integrable classical spin chain, J. Phys. Soc. Jpn. 51 
(1982) 3417-3418. 

[36] V. S. Gerdjikov, M. I. Ivanov and Y. S. Vaklev: Gauge transforma- 
tions and generating operators for the discrete Zakharov-Shabat system, 
Inverse Probl. 2 (1986) 413-432. 

[37] N. Papanicolaou: Complete integrability for a discrete Heisenberg chain, 
J. Phys. A: Math. Gen. 20 (1987) 3637-3652. 

[38] G. M. Pritula and V. E. Vekslerchik: Toda-Heisenberg chain: interacting 
a -fields in two dimensions, J. Nonlinear Math. Phys. 18 (2011) 443-459. 

[39] V. E. Adler and A. B. Shabat: On the one class of hyperbolic systems, 
SIGMA 2 (2006) 093. 

[40] M. J. Ablowitz, B. Prinari and A. D. Trubatch: Discrete and Con- 
tinuous Nonlinear Schrddinger Systems (Cambridge University Press, 
Cambridge, 2004). 

[41] A. Dimakis and F. Muller-Hoissen: Solutions of matrix NLS systems 
and their discretizations: a unified treatment, Inverse Probl. 26 (2010) 
095007. 

[42] G. Baxter: Polynomials defined by a difference system, Bull. Amer. 
Math. Soc. 66 (1960) 187-190. 



65 



G. Baxter: A convergence equivalence related to polynomials orthogonal 
on the unit circle, Trans. Amer. Math. Soc. 99 (1961) 471-487. 

0. O. Vakhnenko: Inverse scattering transform for the nonlinear 
Schrddinger system on a zigzag-runged ladder lattice, J. Math. Phys. 
51 (2010) 103518. 

M. Wadati and T. Kamijo: On the extension of inverse scattering 
method, Prog. Theor. Phys. 52 (1974) 397-414. 

E. Olmedilla: Inverse scattering transform for general matrix 
Schrddinger operators and the related symplectic structure, Inverse 
Probl. 1 (1985) 219-236. 

B. Carl and C. Schiebold: A direct approach to the study of soliton 



equations (2000), see http://apachepersonal.miun.se/~corsch/ 



C. Schiebold: Integrable systems and operator equations, Habilitation 
thesis (2004), see |http: / / apachepersonal.miun.se/~corsch/] . 

T. Aktosun, F. Demontis and C. van der Mee: Exact solutions to the 
focusing nonlinear Schrddinger equation, Inverse Probl. 23 (2007) 2171- 
2195. 

T. Aktosun, T. Busse, F. Demontis and C. van der Mee: Symmetries 
for exact solutions to the nonlinear Schrddinger equation, J. Phys. A: 
Math. Theor. 43 (2010) 025202. 

A. Dimakis and F. Muller-Hoissen: Bidijferential calculus approach to 
AKNS hierarchies and their solutions, SIGMA 6 (2010) 055. 

A. Dimakis, N. Kanning and F. Muller-Hoissen: Bidifferential calculus, 
matrix SIT and sine- Gordon equations, Acta Polytechnica 51 (2011) 
33-37. 

C. S. Gardner, J. M. Greene, M. D. Kruskal and R. M. Miura: Korteweg- 
de Vries equation and generalizations. VI. Methods for exact solution, 
Commun. Pure Appl. Math. 27 (1974) 97-133. 

M. Wadati: On the exact solution of the Korteweg-de Vries equation, 
Science of Light 25 (1976) 37-40. 

T. Tsuchida: N -soliton collision in the Manakov model, Prog. Theor. 
Phys. Ill (2004) 151-182. 



66 



[56] R. Hirota: Exact N-soliton solution of nonlinear lumped self-dual net- 
work equations, J. Phys. Soc. Jpn. 35 (1973) 289-294. 

[57] V. E. Adler and V. V. Postnikov: On vector analogs of the modified 
Volterra lattice, J. Phys. A: Math. Theor. 41 (2008) 455203. 

[58] S. I. Svinolupov: Jordan algebras and integrable systems, Funct. Anal. 
Appl. 27 (1993) 257-265. 

[59] H. Eichenherr and K. Pohlmeyer: Lax pairs for certain generalizations 
of the sine-Gordon equation, Phys. Lett. B 89 (1979) 76-78. 

[60] P. P. Kulish and E. K. Sklyanin: O(N) -invariant nonlinear Schrddinger 
equation — A new completely integrable system, Phys. Lett. A 84 (1981) 
349-352. 

[61] X. Yang and R. Schmid: Backlund transformations induced by symme- 
tries, Phys. Lett. A 195 (1994) 63-73. 

[62] T. Tsuchida and M. Wadati: The coupled modified Korteweg-de Vries 
equations, J. Phys. Soc. Jpn. 67 (1998) 1175-1187. 



67 



